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Dihadron correlations at high transverse momentum in d + Au collisions at y/stm = 200 GeV 
at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy Ion Collider 
(RHIC) . From these correlations we extract several structural characteristics of jets; the root-mean- 
squared transverse momentum of fragmenting hadrons with respect to the jet U (j^), the mean 
sine-squared of the azimuthal angle between the jet axes (sin ((f) jj)}, and the number of particles 
produced within the dijet that are associated with a high-pr particle (dN/dxE distributions). We 
observe that the fragmentation characteristics of jets in d + Au collisions are very similar to those 
in p + p collisions and that there is also little dependence on the centrality of the d + Au collision. 
This is consistent with the nuclear medium having little influence on the fragmentation process. 
Furthermore, there is no statistically significant increase in the value of (sin 2 ((f)jj)} from p + p to d 
+ Au collisions. This constrains the effect of multiple scattering that partons undergo in the cold 
nuclear medium before and after a hard-collision. 
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PACS numbers: 25.75.-q,13.87.-a,24.85.+p 



I. INTRODUCTION 

Jet production in high energy collisions is a useful tool 
to study the passage of scattered partons through a nu- 
clear medium. A dominant hard-scattering process is two 
partons scattering to produce two high-px partons which 
then fragment to produce a dijet. In a nuclear environ- 
ment the partons that participate in the collision can un- 
dergo multiple scattering within the nucleus, potentially 
changing the structure of the dijet. Such changes can 
provide information on the interaction of colored partons 
with the cold nuclear medium. 

Some information on this interaction is already avail- 
able at RHIC energies via the Cronin enhancement of 
the pt spectra pj. In d + Au collisions at RHIC 

Q the cross-section for high-py particle production 
in d + Au collisions is enhanced compared to p + p col- 
lisions, consistent with multiple scattering in the cold- 
nuclear medium increasing the transverse momentum of 
the partons. In this paper we report on a complimentary 
observable to the Cronin effect: the broadening of dijet 
distributions. Such broadening is directly related to the 
additional transverse momentum delivered to the partons 
during multiple scattering, and hence provides a comple- 
mentary handle for comparison between experiment and 
theory. 

Interpreting both d + Au and Au + Au collisions re- 
quires solid knowledge of baseline p + p collisions, espe- 
cially those dijet events at midrapidity that contain two, 
nearly back-to-back, jets produced from a hard (large 
Q 2 ) parton-parton interaction. Experimentally the jets 
are not exactly back-to-back and the acoplanarity mo- 
mentum vector, kx, was measured in p + p collisions at 
ISR energies to have a magnitude, kx, on the order of 

1 GeV/c 5]. This was much larger than the expectation 
that kx was due to intrinsic parton transverse momen- 
tum governed by the hadron size, which would lead to 
kx ~ 300 MeV/c. It was realized early that additional 
gluon radiation either before or after the hard scattering 
will increase the value of kx and the dijet acoplanarity. 

In collisions involving nuclei, multiple scattering within 
the nucleus increases the parton transverse momentum. 
E557 0, E609 0, and E683 9] have all measured an in- 
crease in the dijet acoplanarity with atomic mass of the 
target. In the case of E683, they measured an As depen- 
dence of (k T ) for both 7 + A and it + A collisions. This 
dependence is expected since the number of scatterings 
should be proportional to the length traversed in the nu- 
cleus (L r~j A3), For large A the extracted (k T ) values 
are about 50% above that for collisions with the hydro- 
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gen target, implying that the multiple-scattering effects 
are as important to the broadening of the dijets as are 
the initial state effects at that energy. In the case of p+ A 
reactions, the measured (k T ) values increase more slowly 
than Ai @- Since the (k T ) values show a strong en- 
ergy dependence we need to establish the initial and 
multiple-scattering contributions to (k T ) for p+A reac- 
tions at RHIC energies. The (k T ) values are also known 
to be dependent on the Q 2 of t he p arton-parton interac- 
tion, increasing with rising Q 2 \wi [llj . 

No model is currently available that can reproduce all 
data on the Cronin effect and dijet broadening although 
most include multiple scattering as the underlying mech- 
anism. A recent review 0] considered two large classes 
of models, 1) soft or Glauber scattering where the multi- 
ple scattering is either at the hadronic or partonic level 
and 2) semi-hard multiple scattering where the multiple 
scattering is at the partonic level. 

In both the soft and hard scattering models, the in- 
crease A(k T ) = (k T ) p+ A — {k"x)p+p is proportional to the 
product of the scattering cross section and the nuclear 
thickness function, 

A(k 2 T ) oc v{b, s/a) - 1 = v M s(V^)T A (b) (1) 

where v(b, y/s) is the number of interactions, b is the 
impact parameter of the collision, cms is the multiple 
scattering cross-section, and 7^(6) is the nuclear thick- 
ness function. For the soft scattering models, cms is 
defined to be (Tnn(Vs), the nucleon-nucleon scattering 
cross section, while for the semi-hard models cms is 
a^ H (^/s), the parton-nucleon semi-hard cross-section. In 
the specific case of hard sphere nucleon scattering 0] 

v(b,«Js) = (tnn(Vs)2M 5 "\/ 1 ~ wliere R is the nu ~ 
clear radius, which gives an Ai increase in A(k T ). 

Both types of these models give the same trend in cen- 
trality and the same dependence on the target's atomic 
mass. The difference between them is in the strength of 
the increase with respect to Ta(6) and how this changes 
with beam energy. We will compare the data in this pa- 
per to two specific implementations of the hard-scattering 
models from Qiu and Vitev |14|. and Barnafoldi et 
al. 0. 

An alternative view of the Cronin effect was recently 
proposed by Hwa and Yang . These authors calculate 
the recombination of hard partons with soft partons re- 
leased during the multiple collisions. Because this model 
reproduces the measured Cronin effect at RHIC without 
imparting successive transverse momentum kicks to the 
scattered partons, the authors suggest that there may be 
little to no increase in kx iromp+p to d + Au collisions. 

We also use jet-fragmentation observables to probe 
multiple scattering in cold nuclei, in particular yj (j T ), 
the RMS of the mean transverse momentum of hadrons 
with respect to the fragmenting parton, and the frag- 
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mentation function of the parton, D(z,Q 2 ), where z is 
the fraction of the parton's momentum that a hadron 
carries. If the parton suffers semi-hard inelastic colli- 
sions within the nuclear environment, the parton will 
lose energy and its subsequent hadronization will pro- 
duce fewer high-z fragments and more low-z fragments. 
We cannot directly measure fragmentation functions via 
dihadron correlations, but we do measure the distribu- 
tion of hadrons produced in association with a high-py 
trigger particle. We plot these distributions as a function 
of xe, where xe is defined as 

PT.trig ' PT, assoc , n s 

xe = — pr k (2) 

\PT,tvie\ 

The motivation for the variable xe can be most eas- 
ily seen in the simple case where (z) = 1 for the trig- 
ger particle and the two hadrons are emitted back-to- 
back. In this case, pr,trig is the transverse momen- 
tum of the scattered parton (gT,parton) and for the far- 
side x E = Zassoc = _PT,a SS oc/<7T,parton. Relaxing the as- 
sumption on z t rig, there is still a simple relation be- 
tween xe and z for back-to-back jets at high-p^, where 
xe — Zassoc/ztrig- Hence the dN/dxE distribution for 
hadrons emitted back-to-back from the trigger hadron 
can be related to the fragmentation function: for more 
details see the end of Section [H] 

There is considerable information on xe distributions 
fromp+p collisions. The CCHK collaboration [5J demon- 
strated that the xe distribution scaled, i.e. the distribu- 
tion was approximately independent of pr.trig- Sc aling 
at higher pr,trig was also established by Fisk et al. [T3 , 
and the CCOR collaboration 0,0, providing support 
for the idea that fragmentation of high-p^ partons is in- 
dependent of the momentum of the parton. 

This scaling is however approximate and scaling viola- 
tion was understood by Feynman et al. [f| to be caused 
by the radiation of semi-hard gluons. Scaling violation 
of the fragmentation function D(z, Q 2 ) is now well es- 
tablished experimentally ([20| and references therein). 
For the Q 2 range considered in this paper (10 < Q 2 < 
1000 GeV/c 2 ), the fragmentation functions used in NLO 
calculations [21j drop by 25% for z = 0.6 over the range 
10 < Q 2 < 100 GeV/c 2 . At higher Q 2 the fragmentation 
functions are less dependent on Q 2 , e.g. the fragmenta- 
tion drops by less than 20% at z — 0.6 over the much 
larger range of 100 < Q 2 < 1500 GeV/c 2 . 

In this paper we quantify the extent to which our mea- 
sured x e distributions in d + Au collision scale and com- 
pare the x e distributions to those from p + p collisions 
at RHIC. The goal is to establish whether inelastic scat- 
tering in the cold medium or the recombination mecha- 
nism changes the effective fragmentation function. The 
xe distributions provide a stringent test the recombina- 
tion model from Hwa and Wang j2^|. This model re- 
produces the Cronin effect in d + Au collisions through 
shower thermal recombination mechanism and predicts 
an increase in jet associated multiplicity |22l l23j|. i.e. an 
increase in the near-angle dN/dxE, in d + Au relative to 



p + p collisions. 

The measured xe distributions in d + Au also serve 
as a critical baseline for Au + Au collisions, where the 
strong energy-loss in the dense, hot medium is expected 
to dramatically change the shape of these distributions. 

Our three goals for this paper are 1) to report the 
characteristics of jet structures in d + Au collisions at 
RHIC energies, 2) to establish the extent to which mul- 
tiple scattering changes these structures as a function of 
centrality and by comparison with data from p + p col- 
lisions, and 3) to establish the baseline for jet-structure 
measurements in heavy ion reactions. Any difference be- 
tween jet properties in Au + Au and d + Au collisions 
should be attributable to the hot, dense nuclear matter 
created in the heavy ion collisions. The main results in 
this paper are presented in Section llVI which details the 
measured values of (jy), (sin 2 (0jj)), and the pt, and 
xe distributions from d + Au collisions at ^/snn=200 
GeV. These results are derived from the fitted widths 
and yields of two-particle azimuthal correlations which 
are reported in Section IIVI The experimental methods 
used to obtain these correlations are described in Section 
lllll and the jet quantities we use throughout the paper 
are fully defined in Section ITU 



II. JET ANGULAR (AZIMUTHAL) 
CORRELATIONS 

A. Two particle correlation 

The defining characteristic of a jet is the production 
of a large number of particles clustered in a cone in the 
direction of the fragmenting parton. Traditionally, ener- 
getic jets are identified directly using standard jet recon- 
struction algorithms I* 1 heavy-ion collisions, due 
to the large amount of soft background, direct jet recon- 
struction is difficult. Even in p + A or p+p collisions, the 
range of energy accessible to direct jet reconstruction is 
probably limited to px > 5 — 10 GeV/c, below which the 
jet cone becomes too broad and contamination from the 
'underlying event' background is significant. Jet identi- 
fication is even more complicated for detectors with lim- 
ited acceptance, such as the PHENIX central arms, due 
to leakage of the jet cone outside the acceptance. 

The two-particle azimuthal angle correlation technique 
provides an alternative way to access the properties of 
jets. It is based on the fact that the fragments are 
strongly correlated in azimuth ((f)) and pseudo-rapidity 
(ji). Thus, the jet signal manifests itself as a narrow peak 
in A(j) and Arj space. Jet properties can be extracted on 
a statistical basis by accumulating many events to build 
a A(j) distribution or a A(j) correlation function. Further- 
more, we assume that the soft background is isotropic. 
Distributions in A<fi were initially used in the 70's to 
search for jet signals in p + p collisions at CERN's ISR 
facility 0, Ell [SI • More recently, A<fi distributions and 
correlation functions have been exploited for analysis of 
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FIG. 1: Near-side or single jet fragmentation kinematics in the 
azimuthal plane perpendicular to the beam axis. The angles 
4>tj and 4>aj are the angles (in the plane transverse to the 
beam axis) between the trigger-jet and associated-jet axes, 
respectively. Also shown are the pr vectors for the trigger 
and associated particles, as well as p ut and jr y of the trigger 
particle. 
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FIG. 2: Far-side jet fragmentation kinematics in the az- 
imuthal plane perpendicular to the beam axis. Similar to 
Fig0 the <f> indices aj, jj, and tj denote the angles between 
the associated-jet, jet-jet, and trigger-jet axes. 



jet correlations at RHIC [UHl El El EH • A detailed 
discussion of the two-particle correlation method can be 
found in . These approaches overcome problems due 
to background and limited acceptance, and extend the 
study of jet observables to lower px- 

In the correlation method, two classes of particles are 
correlated with each other - trigger particles, and asso- 
ciated particles. Although the distinction between these 
two classes is artificial, trigger particles are typically se- 
lected from a higher pt range. In this work, we distin- 
guish between two primary categories of correlations: 

Assorted-pr correlation (AC) , where the pr ranges 
of trigger and associated particle do not overlap. 

Fixed-pT correlation (FC) , where the pt ranges for 
trigger and associated particles are identical. 

In this paper, correlations are further categorized via 
a scheme which uses the identity of the trigger and asso- 
ciated particles. Four different types of such correlations 
are presented. Denoting the trigger particle px as PT.trig 
and associated particle pr as pr,assoc , 

• h^-h^ FC 

The pr range of both the trigger and associated 
particle is 1 < pr < 5 GeV/c. 

• h± - ft ± AC 

There are three different selections, 

1. 2.5 < pr.trig < 4 GeV/c with 0.5 < pr.assoc < 
2.5 GeV/c 

2. 4 < p T , trig < 6 GeV/c with 0.5 < p T ,assoc < 4 
GeV/c 

3. 3 < pr.trig < 5 GeV/c with 0.5 < pr.assoc < 3 
GeV/c 

• TT°-h ± AC 

The trigger particle is a neutral pion and the as- 
sociated particle is a charged hadron, where 5 < 
Fr.trig < 10 GeV/c with 1 < pT,assoc < 5 GeV/c. 



B. Extraction of jr, sin 2 ((frjj) from the Correlation 
Function 

In this section, we discuss the framework for the two 
particle correlation method. Figures ^ and El illustrate 
the relation between the two particles and their parent 
jets, for the case when the parents are the same jet or 
the dijet, respectively. The figures also show the rela- 
tionship between jx, sin 2 (0jj) and the kinematic vari- 
ables describing the trigger and associated particle, j't is 
the component of the particle momentum perpendicular 
to the jet momentum. Its projection into the azimuthal 
plane is jx y ■ The quantity p ont (denoted with N or F 
for near- or far-side, respectively) is the component of 
the associated particle's pr that is perpendicular to the 
trigger particle's pr- The vector sum of fcr,i and fcr,2 
produces the dijet acoplanarity and the azimuthal angle 
between the jet axes is <j)jj. 

The RMS value of jr y can be derived from the cor- 
relation functions. For the single jet fragmentation of 
Fig. 2] if we denote A<j>, <fi t j , and 4> a j as the angles be- 
tween trigger-associated, trigger-jet and associated-jet, 
respectively, then the following relations are true: 

Sin(0y) = =2^-trig (3) 

PT.trig 

/ , \ jTy, assoc _ 

SlT!Xy(p a j j — = Xj" iassoc 

PT. assoc 



PT, assoc 

Acp = 4> t .j + 4> aj . 

Assuming <p t j and 4> a j are statistically independent, we 
have (cross terms average to 0) for the near-side, 

(sin 2 Atp N ) = (sin 2 4> t j cos 2 4> a j) + (sin 2 <f> a j cos 2 <f> t j) (4) 

Substituting the sin and cos terms from Eq. [3] into 
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Eq. we obtain the equation for the RMS value of jx v 



V / (^) = ^out,N)/(l + (4)-2( 



2 ^ ^'.trif 



(5) 



where x h = PT, asS oc/PT,trig- 

In the gaussian approximation for the near-side az- 
imuthal distributions, a simple Taylor expansion con- 
nects p ou t with the jet width, a: 



(Pint) = (PT.assoc Sm2 A <i 



(PT, 
(Pi, 



.) [sin(A</> 2 ) - 



(a</> 4 ). 



(0) 



Since Eq. contains the variable a^.trig that depends 
on j Ty , we should calculate \J(jr y ) heratively. In cases 
when trigger and associated particle px are much larger 
than the typical jx value, the near-side jet width is 
0. Hence Eq. [3] can be simplified as, 

,assoc/ 



small and x 



j.tng 



(7) 



a N 



(PT,trig)(p 



'T.assoc/ 



\T<PT, 



trig/ 



T,assoc/ 



Since jx is the projection of hadron px perpendicular 
to Pxjct, jr y is necessarily less than px- So, for any 



given px range, there is always a upper kinematic cut off 
on the jx y distribution. This effect, known as the seagull 

effect, leads to a reduction on the observed «/ 1 (Jt v ) f rom 
the expected value. It is important at low PT,trig and 
becomes negligible once pr.trig ^> \J(jx v )- The seagull 
effect can be parameterized and removed from the jx y 
values |3^|. 

For the far-side correlation from Fig. [2] we have 



7T — A<j}p = 

sin(A0i?) = 



<Ptj + 4>aj 
Pout,F 

PT,assoc 



+ ■ 



(8) 



where ^jj is the azimuthal angle between the two jet 
axes. Expanding sin 2 A<f)p and dropping all cross terms 
(which average to 0), we get 



(sin 2 A4>p) = ^(sin 4> t j cos 4> a j cos 4>jj) 2 ^ + 
(sin 4> a j cos 4>tj cos cfijj )' 



(sin cf>jj cos 4> a j cos 4> tJ )' 
(sin (j> tj sin cj> aj sin <j>jj) 



(9) 



We substitute Eq. 0]to get 



(sin 2 A</>j?) = (sin 2 A(f> N ) (cos 2 + (cos 2 A(fj N ) (sin 2 0^) 

= (sin 2 A4>n) (l — sin 2 0^-) + (l — sin 2 A^jy) (sin 2 0^-) 

= (sin 2 A0iv) + (sin 2 <^-) - 2 (sin 2 A^jv) (sin 2 (10) 



Collecting terms in <pjj produces 



sin 



) = 



(sin 2 A(f> F ) - (sin 2 A<j> N ) 
1 - 2 (sin 2 A(j) N ) 



(11) 



Note that since is the azimuthal angle between the 
jet axes, sin 2 (jfijj) is one measure of the extent to which 
the jets are not back-to-back, and is hence a quantity that 
is sensitive to any additional scattering in d + Au colli- 
sions. We express the right side in terms of the observ- 
ables ct/v and <tf, the rms widths of distribution that we 
measure by expanding the sine term 

(sin 2 A(f>) = a 2 - ct 4 + 2/3ct 6 (12) 

which is good to 2% for rms widths at 0.5 rad and good 



to 0.6% for rms widths of 0.2 rad. Therefore 



sm 



) = 



+ 2/3<4 ) - (a% -a% + 2/3a%) 



1 - 2 « - a% + 2/3a 



(13) 

The right hand side is now in terms of experimental ob- 
servables which we will use to extract sui((j)jj). 

We have attempted to extract kx from sin 2 (<fijj ) . This 
requires assumptions on the scattered quark distribution, 
the magnitude of the momentum asymmetry between the 
partons due to the kt-kick, as well as the detailed shape 
of the fragmentation function. The current paper is fo- 
cussed on the comparison between p+p and d + Au colli- 
sions, which can be made with sin 2 (<f> jj). Hence we leave 
the extraction of kx to future work. 

In this paper we report the RMS values of jx and 

sin 2 '{(f) jj): where y/ (j T ) = \pl \J{jx y ) l n tne literature, a 
jx value is sometimes reported as the geometrical mean, 
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(|j T J)- The relation to the RMS value is (Ij'tJ) 



C. Conditional Yields 

We also present in this paper the associated yield per 
trigger particle, referred to as the conditional yield (CY) , 
as a function of pr and xe- The CY is the number of 
particles produced in the same or opposite jet associated 
with a trigger particle 



PHENIX Detector 

PC.3 



CY( PT ) 
CY(x E ) 



1 dN h 

1 dN h 
Ntrig dx E 



(14) 
(15) 



and can be directly extracted from the measured Gaus- 
sian yields in the correlation functions. 

To emphasize the importance of the CY, we note that 
it is related to the single- and two-particle cross sections: 



CY = 



dpadpb 



' da 
dpr 



(16) 



The interpretation for the two-particle cross section de- 
pends on whether one is studying the near- or far-side 
jet correlations. The conditional yield for particles from 
the near-side jet depends on the dihadron fragmentation 
function, while the conditional yield from the far-side 
jet depends on two independent fragmentation functions: 
one parton fragments to produce a hadron with pT,trig ; 
while the other scattered parton on the far-side fragments 
to produce a hadron with pr.assoc- For the far-side con- 
ditional yield at high-p^, %e — -Zassoc/^trig (see Eq. 
Hence, d(xs) — d(z assoc ) / z tr i g and the slope of the far- 
side CY(xe) is ^trig times the slope of the fragmentation 
function D(z). 



III. EXPERIMENT AND DATA ANALYSIS 
A. Data Collection 

The data presented in this paper were collected by the 
PHENIX experiment at the Relativistic Heavy Ion Col- 
lider during the d + Au and p + p run of January - 
May 2003. During that time integrated luminosities were 
recorded of 2.7 nb _1 for d + Au collisions and 0.35 pb _1 
for p + p collisions each at ^/snn=200 GeV. 

The PHENIX detector consists of two central spec- 
trometer arms, two forward muon arms and several global 
detectors used for triggering, vertex detection, and cen- 
trality selection. This analysis utilizes the two central 
spectrometer arms that each cover a region of |?7| < 0.35 
units of pseudorapidity and 90° in azimuth. The spec- 
trometer arms are not exactly back-to-back in azimuth, 
so while there is large acceptance for the detection of two 




West 



Beam View 



East 



FIG. 3: (Color online) The two central spectrometer arms of 
the PHENIX experiment used to collect the charged hadron 
and charged and neutral pion tracks. 



particles separated by 180° there is also finite acceptance 
for two particles separated by 90°. Figure|3|shows a beam 
cross-section view of the PHENIX central spectrometer 
arms. A complete overview of the whole PHENIX detec- 
tor is found in reference [sH ■ In this section we will only 
focus on those subsystems relevant for the analysis of the 
dihadron data. 



1. Global Event Characteristics 

For event characterization the Beam-Beam Counters 
(BBCs) |13 are utilized. The BBCs are sets of 64 
Cherenkov counters placed symmetrically along the beam 
line, covering 3 < |^| < 3.9 units of pseudorapidity and 
located 144 cm from the center of the interaction region. 
The BBCs determine the event vertex and the initial col- 
lision time, to, from the time difference between parti- 
cles reaching each BBC. For this analysis we include only 
events with an offline cut of |z ve rtex| < 30 cm. 

The BBC facing the direction of the Au beam was 
used to determine the centrality. Fig. 0] shows the BBC 
charge distribution and the centrality classes used in this 
analysis. The centrality is defined as 

%Centrality = 88.5%(1 - frac(Q BB c)) (17) 

where frac(QeBc) is the fraction of the total BBC charge 
distribution integrated from zero to Qbbc and 88.5% is 
the efficiency of the minimum bias trigger. This central- 
ity can be related to the mean number of Au participants, 
(-^Vpart), and mean number of collisions, (N co \\). To de- 
termine the centrality we model the BBC charge distri- 
bution as a negative binomial distribution with a width 
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TABLE I: Table of the mean number of collisions for d + 
Au , iVcoii, versus the percentage of the total inelastic cross 
section and the nuclear overlap function Ta(o). 



Percent <7j ne i 


(iVcoll) 


(T^Xmb- 1 ) 


0-20% 


15.4 ± 1.0 


0.367 ±0.024 


20-40% 


10.6 ±0.7 


0.252 ±0.017 


40-88% 


4.7 ±0.3 


0.112 ±0.007 



and mean proportional to -/V part . So, for a given central- 
ity, there are several negative binomial distributions (de- 
fined by A^part) that contribute to the overall distribution 
and as such N paTt is not uniquely defined. We calculate 
a weighted average of N pait , where the weight is given by 
the negative binomial distribution for a given -/V part and 
the probability for having a collision with N pa _ vt . The lat- 
ter probabilities were computed using a Glauber model, 
with a Hulthen wave function for the deuteron and an 
inelastic cross section of 42 mb. Finally, the (N co a) was 
determined for a given (N pax t) from the same Glauber 
model. The resulting centrality bins and (N co \i) used in 
this analysis are outlined in Table [I] 

The dihadron events were recorded using several dif- 
ferent Level- 1 triggers. The minimum bias trigger re- 
quired at least one hit in each of the BBCs and that 
the collision vertex (computed online) satisfies |z V crtox| < 
75 cm. It was sensitive to 88.5% of the inelastic d 
± Au cross section. PHENIX also employed a se- 
ries of Level-1 triggers to select electrons, photons and, 
with lower efficiency, high-p-r hadrons. These trig- 
gers are called the ERT triggers and they utilized the 
Ring Imaging Cerenkov (RICH) for electron identifi- 
cation, together with the Electromagnetic Calorimeter 
(EMC) j3y|, which consists of 8 sectors, 6 of which are 
Lead-Scintillator (PbSc) sampling calorimeters and 2 are 
Lead-Glass (PbGl) Cherenkov counters. The EMC has 
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FIG. 4: (Color online) Total charge distribution on the Au- 
going side Beam-Beam Counter (BBC) for d ± Au collisions 
and the centrality selection (see Table HJ. 



excellent timing and energy resolution for electromag- 
netic showers. The ERT triggers were produced by sum- 
ming signals from tiles, where a tile was 4x5 photo- 
multipliers (PMTs) in the RICH and either 2x2 or 4x4 
PMTs in the EMC. 

The electron trigger was defined by the coincidence be- 
tween the minimum bias trigger and the RICH and EMC 
2x2 trigger where the threshold for the RICH tile was 3 
photo-electrons and the EMC threshold varied between 
400-800 MeV. Three different thresholds were available 
for the 4x4 photon triggers. These thresholds differed 
between the PbGl and PbSc and varied within and be- 
tween the p ± p and d ± Au runs. The lowest threshold 
setting (1.4 GeV-2.8 GeV) was most sensitive to hadron 
showers in the EMC. The threshold values and rejection 
factors (rejection = N ln i n BiasEvcnts/N t riggcrEvcnts) for the 
ERT triggers, in coincidence with the minimum bias trig- 
ger, are given in Table [H] The h — h correlations 
use only the minimum bias triggered data, while the 
7T° — correlations use only the ERT photon triggers. 
The tt ± — h correlations use the minimum bias, ERT 
photon and ERT electron triggers. A detailed knowl- 
edge of the ERT trigger efficiency is not necessary, since 
we present the conditional yield distributions per trigger, 
for which this efficiency cancels out. 



2. Tracking and Particle Identification 

In this section we discuss the tracking and identifi- 
cation of the particles used in the different correlation 
analyses. There are three types of particles included: 
Charged hadrons are used in all analyses, neutral pions 
are used as trigger particles for the ir° — h correlations 
and charged pions are used as trigger particles for the 
7T^~ — hF 1 correlations. 

Charged hadron tracks are measured outside the 
PHENIX central magnetic field by the Drift Chamber 
(DC), located 2.0 m from the vertex, and two layers 
of multi-wire proportional chamber (PCI and PC3), lo- 
cated 2.5 and 5.0 m, respectively, from the vertex |37| . 
The DC determines the momentum and the azimuthal 
position of the track, while PCI determines the polar 
angle |38j. The momentum resolution is determined to be 
O.7%01.1%p (GeV/c)@. Tracks are confirmed by re- 
quiring that an associated hit in PC3 lies within a 2.5 a 
(for h —h ) or 3<r (for TT — h ± ) matching window in both 
the 4> an d z directions. This cut reduces the background 
from particles not originating in the direction of the ver- 
tex. The remaining background tracks are mainly decays 
and conversion particles 39] . The background level for 
single-tracks is less than 5% below 3 GeV/c, increasing 
to about 30% at 5 GeV/c. However, the background is 
smaller for high-px triggered events (see Section IIII B|l . 
The charged particle tracking efficiency for the active re- 
gion of the DC, PCI and PC3 is better than 98%. Since 
we perform a pair analysis, the two track resolution is 
important. For the DC, the two track separation is bet- 
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TABLE II: EMC threshold and rejection factors for the electron and photon ERT triggers in coincidence with the minimum 
bias trigger for p+p and d + Au. The photon triggers are defined by the energy sum of 4x4 PMTs in the EMC above threshold. 
The electron trigger requires the coincidence of the RICH trigger (threshold of 3 photo-electrons for both p+p and d + Au runs) 
and the energy sum of 2x2 PMTs in the EMC above threshold. 





p- 


~P 


d + A 


u 


P + P 


d + Au 




PbSc threshold 


PbGl threshold 


PbSc threshold 


PbGl threshold 


rejection 


rejection 


Gammal 


2.1 GeV 


2.1 GeV 


2.8 GeV 


3.5 GeV 


400-1200 


125-300 


Gamma2 


2.8 GeV 


2.8 GeV 


2.8-3.5 GeV 


3.5-4.2 GeV 


1500 - 3100 


450-900 


Gamma3 


1.4 GeV 


1.4 GeV 


2.1 GeV 


2.8 GeV 


70-160 


15-60 


Electron 


0.4-0.8 GeV 


0.4-0.8 GeV 


0.6-0.8 GeV 


0.6-0.8 GeV 


5-1200 


30-170 



ter than 1.5 mm, while at PCI it is 4 cm and at PC3 it 
is 8 cm. 

Neutral pions are detected by the statistical recon- 
struction of their 77 decay channel. These decay pho- 
tons are detected by the EMC and identified by their 
time-of-flight (TOF) and shower shape. The electromag- 
netic shower shape is typically characterized by the x 2 
variable 13a. 



x1(T 



(Ef 



E 



pred\2 



(18) 



where £™ eas is the energy measured at tower i and Ef led 
is the predicted energy for an electromagnetic particle of 
total energy X^-^ mcas - This \ 2 value is useful for the 
discrimination of electromagnetic from hadron showers. 
The x 2 an d TOF cuts used give a very clean sample of 
photons with contamination of other particles at <C 1%. 

Using pairs of photons that pass these EMC cuts, we 
create the invariant mass spectra for each photon pair 
Pt- A sample invariant mass distribution with a S/B of 
approximately twelve is given in Fig. [5] The background 
distribution can be reproduced by mixing clusters from 
different events and normalizing that distribution to the 
real event distribution outside the ir° mass region. The 
peak position and width of the invariant mass distribu- 
tion were parameterized as a function of pair px, in order 
to select 7T° candidates from a region of invariant mass 
within 2er of the peak position. The S/B for a tt° with 
Pt > 5 GeV is 10 — 20, increasing as a function of pr- 
There is a slight dependence on centrality with the tt° 
S/B decreasing with increasing centrality. 

PHENIX identifies high momentum charged pions 
with the RICH and EMC detectors. Charged parti- 
cles with velocities above the Cherenkov threshold of 
7th = 35 (CO2 radiator) emit Cherenkov photons, which 
are detected by photo-multiplier tubes (PMTs) in the 
RICH 40]. This threshold corresponds to 18 MeV/c 
for electrons, 3.5 GeV/c for muons and 4.9 GeV/c for 
charged pions. In a previous PHENIX publication [3^], 
we have shown that charged particles with reconstructed 
Pt above 4.9 GeV/c, which have an associated hit in 
the RICH, are dominantly charged pions and background 
electrons from photon conversions. The efficiency for de- 



1500 



> 

a) 

O 1000 

S * 

5 500 

Z 

T3 

a 



d+Au 0-20% central 
6 < p T < 7 GeV/c 

S/B = 12.5 



0.1 0.2 0.3 

[GeV/c 2 ] 



0.4 



FIG. 5: Measured 77 invariant mass distribution for 6 < pr < 
7 GeV/c in central d + Au collisions. The peak is fitted with a 
gaussian to extract the centroid mass and a. The S/B within 
2 a of the centroid ranges from ~ 6 at 3 GeV/c up to 20 at 
~ 8 GeV/c. 



tecting charged pions rises quickly past 4.9 GeV/c, reach- 
ing an efficiency of > 90% at pt > 6 GeV/c. 

To reject the conversion backgrounds in the pion candi- 
dates, the shower information at the EMC is used. Since 
most of the background electrons are genuine low pt par- 
ticles that were mis-reconstructed as high pt particles, 
simply requiring a large deposit of shower energy in the 
EMC is very effective in suppressing the electron back- 
ground. In this analysis a momentum dependent energy 
cut at EMC is applied 



E > 0.3 + 0.15p T 



(19) 



In addition to this energy cut, the shower shape infor- 
mation |36j is used to further separate the broad hadronic 
showers from the narrow electromagnetic showers and 
hence reduce the conversion backgrounds. In this analy- 
sis we use the probability (prob) calculated from the % 2 
value (Ea. ll8|) for an EM shower. The probability values 
range from to 1 , with a flat distribution expected for an 
EM shower and a peak around for a hadronic shower. 
Fig.[S]shows the probability distribution for the pion can- 
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didates and electrons, normalized by the integral, where 
the pions candidates were required to pass the energy 
cut and the electrons were selected using particle ID cuts 
similar to that used in Ref. . Indeed, the electron dis- 
tribution is relatively flat, while the charged pions peak 
at 0. A cut of prob < 0.2 selects pions above the energy 
cut with an efficiency of > 80%. Detailed knowledge of 
the pion efficiency is not necessary, since we present in 
this paper the per trigger pion conditional-yield distribu- 
tions, for which this efficiency cancels out. 
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FIG. 6: (Color online) The probability distribution for 
charged pion candidates and electrons. The histogram in- 
tegral has been normalized to unity. 



Since the energy and prob cuts are independent of each 
other, we can fix one cut and then vary the second to 
check the remaining background level from conversions. 
The energy cut in Eq. ^| is chosen such that the raw 
pion yield is found to be insensitive to the variation in 
prob. Fig.|7|shows the raw pion spectra for ERT triggered 
events as function of px, with the above cuts applied. 
The pion turn on from 4.9 — 7 GeV/c is clearly visible. 
Below pt of 5 GeV/c, the remaining background comes 
mainly from the random association of charged particles 
with hits in the RICH detector. The background level is 
less than 5% from 5 — 16 GeV/c, which is the px range 
for the charged pion data presented in this paper. 



B. Data Analysis 

In this section we outline the method to obtain corre- 
lation functions and distributions. From these we extract 
the jet shapes and yields outlined in Section [HI For the 
extraction of the jet-yield from the azimuthal distribu- 
tions we discuss how we obtain the absolute normaliza- 
tion of the distribution, while for the jet shape proper- 
ties, jx v and (sin 2 <pjj), the absolute normalization is not 
necessary. 
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FIG. 7: (Color online) The raw charged pion transverse mo- 
mentum spectrum, with the final cuts applied. The level of 
the remaining background is estimated from an extrapolation 
from low-pr and is shown as a black line. 



Correlation Functions 



Azimuthal correlations functions are generally defined 



as 



C (A<£) cx 



jVcor (A0) 

AW (A<j>) 



(20) 



Similarly, one can also define the correlation function in 
pseudo-rapidity, 



C7(Ar ? )cx A "" 



(A V ) 



A mix (At?) 



(21) 



The same-event pair distribution, N COT (A<f>) or A cor (A^), 
is constructed for trigger-associated particle pairs. The 
mixed-event pair distribution, A m i X (A</>) or A m i X (A?y), 
is determined by combining trigger particles with associ- 
ated particles from randomly selected events. 

This definition of the correlation function relies on the 
fact that detector acceptance and efficiency cancels. It is 
therefore important that the pair efficiencies of the aver- 
age mixed event background and the average foreground 
distributions are the same. For this reason we gener- 
ate mixed event distributions only for events with similar 
centrality and event vertex. More precisely, mixed events 
were required to match within ± 10% centrality and the 
event vertices were also required to be within ± 3 cm. 
For h — h ± correlations the real and mixed events are 
minimum bias data. For 7r° — h correlations the real 
and mixed events are ERT-triggered data. For tt^ 1 — h 
correlations the real events are ERT-triggered and mini- 
mum bias data while the mixed events mix ERT-triggered 
events with minimum bias events. 

For h — h and ^ — h correlations, due to finite 
two track resolution for charged particles at the DC and 
PC, the reconstruction efficiency for same-event charged 
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track pair drops at small Acj> and Ar]. To minimize the 
difference of the pair efficiency between N COI (A<f)) and 
iV m i x (A<^), the pairs are required to have a minimal sep- 
aration of about two times the resolution at the various 
tracking detectors. This corresponds to about 0.28 cm, 
8 cm and 15 cm at the DC, PCI and PC3, respectively. 
However, these pair cuts are not required for ir° — cor- 
relations, because different detector subsystems are used 
for reconstructing trigger- ir° and the associated charged 
tracks as outlined earlier. 

Given the similarity of the analysis techniques between 
A(f> and A-q correlations, in this paper we focus on the A<p 
correlation. The A(f> correlation functions are obtained 
with two different normalizations. For ft, ± — hr 1 assorted 
correlations, the correlation function is area normalized 



C norm (A0) 



N COI (Acj>) J dA0{N n 



<(A0)) 



7V mix (A0) / dA0(iV cor (A0)) 



(22) 



The details concerning this normalization are discussed in 
Section TlV Al The second normalization is used in both 
the 7T° — h and ^ — correlations. It was shown in 
Ref [3^ that the CY can be derived from the measured 
correlation function with an appropriate normalization, 



1 dN _ i? A „ N COT (Act)) 



N° dAcf) N tlig e JUt^dMl 

J dA<f>N Inix (A4 



(23) 



where N® Ti& and iVtrig are the true and detected number 
of triggers respectively, and e is the average single parti- 
cle efficiency for the associated particles in 2n in azimuth 
and ±0.35 in pseudo-rapidity. Ra v corrects for the loss 
of jet pairs outside a given Arj acceptance, determined 
by PHENIX's finite acceptance in 77. This second nor- 
malization is defined so that the integral of the resulting 
correlation function should be No/N® lig , the total num- 
ber of pairs per trigger particle in a given azimuthal and 
eta range. 

For the normalization in Ea. l23l there are two separate 
efficiencies that must be determined, the Arj correction 
and the single particle efficiency. The near-side correla- 
tion has a well-behaved peak around Arj — 0. As we 
show in Section II V Al the near-side jet width in A<fi and 
A77 are consistent with each other within errors. So we 
correct the near-side yield to the full-jet yield assuming 
the shape of the jet is Gaussian and the widths ar e eq ual 
in Acf> and Arj. This correction, according to Ref. |32j | , is 



RAt) — 



1 



J-o.rdAn 



z(A v ) 



(24) 



where acc(Ar)) represent the PHENIX pair acceptance 
function in |A?j|. It can be obtained by convoluting two 
flat distributions in \rj\ < 0.35, so acc(Ar]) has a sim- 
ple triangular shape: acc(Arj) = (0.7— \Arj\) /O.7. The 
PHENIX single particle acceptance is flat in r\ to within 
5%. 

In the far-side the jet signal is much broader than 
the PHENIX acceptance due to the broad range of 
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FIG. 8: The correction factor Ra^ plotted as function of jet 
width for near-side (dashed line) and far-side (solid line). 



momentum-fraction x of the partons that participate in 
the hard-scattering. In fact, we studied the far-side jet 
shape for tt^ 1 — correlation (Fig. 117b ) and found the 
true jet correlation strength is almost constant in the 
PHENIX pair acceptance |A?7| < 0.7. Based on that, 
we assume that the far-side jet strength is constant and 
correct the far-side yield to the corresponding accessible 
pair range of \Ar]\ < 0.7, 



R 



A i) 



2 x 0.7 



0.7 J 



(25) 



R^ equals two, because the pair efficiency has a trian- 
gular shape in |A?7| < 0.7, which results in 50% average 
efficiency when the real jet pair distribution is flat in 
|A?7| < 0.7. Figure |H1 shows the correction factor R& v as 
a function of jet width. The typical range of the near- 
side jet width in all analysis (see Section TlV A(l is below 
0.5 radian. The maximum correction is about a factor of 
2 for the near-side jet. 

The single particle efficiency for associated particles, 
e, includes detector acceptance and reconstruction effi- 
ciency. It is evaluated in a way similar to previously pub- 
lished Au + Au [33 and d + Au 2j analyses. However, 
the jet associated charged hadron spectrum in d + Au is 
much flatter than the inclusive charged hadron spectra 1 , 
so the corrections due to momentum scale and momen- 
tum resolution are much smaller than that for inclusive 
charged hadron. For the same reason, the background 
contamination at high pt, mainly coming from decay 
and photon conversions which are falsely reconstructed 
as high pt tracks 0, H|J, is also reduced. We studied 



1 For example in Fig. 1211 the jet associated yields decreases by a 
factor of 100 from 0.5 to 5 GeV/c. However the typical single 
inclusive hadron spectra decrease by a factor of 100000 2] . 
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both effects using a full GEANT simulation of PYTHIA 
events through PHENIX detectors. The jet associated 
yields were extracted in the same manner as for the real 
data analysis. By comparing it with the input jet as- 
sociated yield spectra, we can quantitatively study the 
effect of momentum smearing and high px background 
contamination. The corrections due to momentum scale 
and resolution are found to be less than 5% with 3% sys- 
tematic errors. For high-py triggered events, the back- 
ground contamination to the associated charged hadrons 
is found to be 5% independent of pt from 1 — 5 GeV/c. 



2. Extracting jet properties 

The normalized correlation functions and conditional 
yield distributions are both fitted with a sum of two 
Gaussians to extract the jet widths and the conditional 
yield of hadrons in the near-side (Ac/) ~ 0) and far-side 
(Ac/> ~ 7r). The fit for the normalized correlation func- 
tions is described in Section IIV Al For the conditional 
yield, we fit with the following function 



1 clN Yield N Yieldr 

= jd -\- c N -\- e F 

Ng ig dAcp V2^<J N V2^ct f 

(26) 

where B reflects the combinatoric background level in 
the real distribution relative to the mixed distribution 
and the other two terms represent the near-side jet and 
far-side jet signal, respectively. The resulting widths, ctjv 
and <7f, are then used to calculate the jet shapes via 
Eas.l7landrHfl For each choice of trigger and associated 
particle pr range, Yield at and Yields directly reflects the 
jet associated yield, dN /dpr, aS soc at the near and far-side, 
respectively. 

Two methods were used to calculate the dN/dxE dis- 



ci 



corre- 



tribution. The first method was used for tt 
lations. Since these correlations are binned in px there 
is a distribution of xe for each trigger-associated px bin. 
This distribution is approximately Gaussian. The fit- 
ted peak value is used as the bin center of the dN/dxE 
distributions and the fitted gaussian width is used as the 
horizontal error bar. To estimate the bin- width in xe, we 
used the definition Eq. [21 which can be written as (ignor- 
ing the sign) x E = PT.assoc cos(A0)/pT )trig . We estimate 
cos(Acj)) ps 1 and write the bin-width as 



Ax, 



^ m ax mm 

^T, assoc ^T, assoc 



(PT.tns) 



(27) 



mm ^max 
assoc T.assocJ 



where we have an associated px bin from [p™ 
and a trigger px bin with a mean (pr,trig)- 

The second method is adopted by it*- — analysis. 
It is statistically based and can be used to calculate the 
distribution for any pair variable pr, trig, PT,assoc, Ac/), Ar], 
x e, Pout etc. In the following we show two examples: 
the dN/dAcj) and dN/dxE distributions. For each pair 
we calculate the Acf> and x e value, then from Eq. [231 we 



calculate the same correction factor that was used for the 
dN/dA</) distribution. 



w(Acp) 



Rat 



1 



7V trig£ 27rjV mix (A0) 



(28) 



If this weight is used to fill the Ac/> histogram for the 
real and mixed distribution, we obtain the CY for the 
true real pairs, 



1 dN 
K~dA^ 



= A(/) roa iw(A0 rcal ) 



(29) 



and for the mixed pair the sum should be one by defini- 
tion, 

background(Acj)) = Ac/> m i VL w(Ac/) m i X ) = 1 . (30) 



Thus the jet signal can be extracted as 
1 dN' et 

N o ri& rfA ° = E A( ^ w ( A< ^)- ( 31 ) 

(Acp mix ) , 

i 

consistent with Eg 1261 

When this weight is used to fill the xe histogram for 
both real and mixed distributions, we obtain the dN/dxE 
by subtracting the mixed xe distribution from the real 
xe distribution, 



1 rl /V jet 



N° dx E 



(32) 



Equation 1311 is rather trivial, because the weighting 
procedure is equivalent to Eq. 1231 for which we know the 
shape of the distribution (Eq. I26|) . But the advantage of 
the weighting procedure is that it allows for the deter- 
mination of the absolute background pair distribution in 
any pair variables. 

Similarly, the statistical method is used to extract the 
PT,assoc and p ou t spectra as 



1 dNt_ = ^p-a! (A ^ cai) 
A t U rig dpT, assoc ' 



(33) 



Bx^£ soc ,(At) 



1 dN^ ct 

= Aureal) - (34 ) 

trig u Pout 



By construction, the integral of the jet yield should be 
conserved independent of the pair variable used, i.e.: 
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dA<j> 



dN 
dA(f> 



dxr 



I 



dp' 



T.assoc 



dN 3 



id 



dp 



'T, assoc 



da;? 



= / dp 



(35) 



id 



dp 



C. Systematic Uncertainties 

The correlation analyses presented here consist of sev- 
eral steps ranging from the generation of correlation func- 
tions to the extraction of the final physics variables (Jt, 
sin((f)jj), per trigger yields, etc.) from these correlation 
functions. Systematic error estimations for each of these 
steps have been evaluated and combined to determine 
the overall error quoted for each measurement. All er- 
rors quoted are maximum extent. 

Systematic errors associated with the generation of 
correlation functions can result from shape distortions 
in either the foreground or the background distributions. 
These distortions can arise if the requisite quality cuts 
(see section UlI All are not stable. In order to minimize 
such errors, the track-pair and quality cuts were assigned 
such that the correlation functions were essentially in- 
sensitive to reasonable cut variations. Systematic er- 
rors associated with such cut variations are estimated 
to be less than 4%. A further source of systematic er- 
rors is related to the efficiency of the background rejec- 
tion when requiring a confirmation hit in the outer pad 
chamber. The yields have been corrected for remaining 
background. The systematic error on the background es- 
timate is ~ 3% for tracks with a transverse momentum 
(p T ) < 4 GeV/c and w 7% for particles with 4 < p T < 6 
GeV/c. For the calculation of the conditional yields, the 
systematic error is dominated by the uncertainties asso- 
ciated with the determination of the efficiency corrected 
single particle yields. These systematic errors have been 
estimated to be ss 10% as obtained from Ref 0. This 
error has two parts; the normalization error includes the 
error on PC3 matching and active area. The momentum 
smearing error includes contributions from momentum 
resolution and momentum scale. 

A separate error is estimated for tt° — correlations 
due to the background contamination of the 7r°s within 
the mass bin. To estimate the width and yield contribu- 
tion of the background 77 pairs, we created correlations 
of 77 outside the tt° mass with hadrons. From these 
we extrapolated the background contribution at the ir° 
mass. These systematic errors are pr-dependent. For the 
near and far angle width the variation is 1 — 3%, the near 
yield variation is 1%, the far yield variation is 1 — 5% and 
increases with increasing px- 

The event mixing technique has been used to correct 
for the limited detector acceptance and inefficiency. In 
addition, the CY has been corrected for limited A77 cov- 
erage. To cross check these procedures we have run a de- 
tailed simulation using the PYTHIA event generator (4^ 



coupled to a single particle acceptance filter that ran- 
domly accepts charged particles according to the detec- 
tor efficiency. In the following, we shall use tt^ 1 — h as 
an example for this cross check. Figure [5] shows a typi- 
cal PHENIX two dimensional single particle acceptance 
used in this analysis. 




4)(rad) 



FIG. 9: (Color online) A typical PHENIX single particle 
acceptance for charged hadrons. 



We generated 1 million PYTHIA events, each required 
to have at least one > 6 GeV/c charged pion. To speed up 
the event generation, a Q 2 cut of 100 GeV 2 on the under- 
lying parton-parton scattering is required. These events 
were filtered through the single particle acceptance fil- 
ter. As an approximation, we ignore the pr dependence 
of acceptance. The same event and mixed pair A(f> dis- 
tributions were then built by combining the accepted ir^ 
and charged hadrons. The jet width and raw yield were 
extracted by fitting the -j^; I ^jff^ with a constant plus 
double gaussian function. The raw yields were then cor- 
rected via Eq.|31to full jet yield for the near-side and the 
yield in |At/| < 0.7 for the far-side. We also extracted the 
true CY and jet width without the acceptance require- 
ment. The comparison of the CY and jet width with and 
without the acceptance requirement are shown in Fig llOl 
The trigger particles are ^ with 6 < pr.trig < 10 GeV/c, 
the associated particles are h ± . In the near-side, the 
corrected yield (top left panel) and width (bottom left 
panel) are compared with those extracted without ac- 
ceptance filter. In the far-side, the yield corrected back 
to |Ary| < 0.7 (top right panel) and the width (bottom 
right panel) are compared with those extracted without 
the acceptance filter. The data requiring the acceptance 
filter are always indicated by the filled circles, while the 
expected yield or width are indicated with open circles. 

The agreement between the two data sets can be better 
seen by plotting the ratios, which are shown in Fig. 1771 
The yields agree within 10% and the widths agree within 
5%. Since yj (jr), (sin 2 (</> J j)) are derived from the jet 
widths, the agreement in width naturally leads to the 
agreement in the \J (j^) and (sin 2 '{4>jj))- O ne notices 
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FIG. 10: (Color online) The comparison of near-side yield 
(top left panel), near-side width (bottom left panel), far-side 
yield (top right panel) and far-side width (bottom right panel) 
as function of pr of charged hadrons. These are obtained 
for 71"* — h correlation from PYTHIA, with a trigger pion 
of 6 — 10 GeV/c. The filled circles represent the quantities 
calculated with PHENIX acceptance filter. 
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FIG. 11: (Color online) The ratio of the jet width or corrected 
yield obtained from the event mixing method to those without 
the acceptance filter. 



that there are some systematic differences in the com- 
parison of the yield at low pr.assoc • This might indicate 
that the Gaussian assumption is not good enough when 
the jet width is wide and the extrapolation for \ Ai]\ > 0.7 
become sizeable. (At px.a SSO c = 0.5 GeV/c, the jet width 
<7at = 0.5 (rad) and the extrapolation is about 20%.). 

The approximations in the formulas used to extract jx y 
and sin 2 (<f)jj ) are used to estimate the systematic error on 
these quantities. We estimate the systematic uncertainty 
in the formulation at the level of 5% for the y/TJj) and 



Table IIIII summarizes the systematic errors for the 

extracted widths, y/TJx) an< ^ \J ( s ^ 1 ° 2 ( ( f > jj)) > while Ta- 
bles lIVllVl and lVII summarize the list of systematic errors 
on the CY for the hadron-hadron, neutral pion-hadron, 
and charged pion-hadron correlations, respectively. Ta- 
ble [^n] outlines the systematic errors on the p out extrac- 
tion from pion-hadron correlations. 



IV. RESULTS 

We present the minimum bias and centrality depen- 
dent results on extracted jet widths and yields in Sec- 
tion II V Al which are used in Section II V Bl to calculate 
quantities describing the jet-structures: the values of 
\J (j T ), (sin 2 (0jj )) , and jet fragmentation conditional 
yields dN/dpx and dN/dxE- The minimum-bias d + 
Au results are compared with results from p + p in Sec- 
tion llV Cl to establish the extent of effects due to medium 
modification in d + Au with as much statistical preci- 
sion as possible. The d + Au centrality dependence of 
the derived quantities is presented in Section ITVDI This 
provides a larger lever-arm in nuclear thickness function, 
at the cost of dividing the available minimum-bias data 
into different centrality bins. 



A. Correlation Functions, Widths and Yields 

The baseline data from which jet structures are ex- 
tracted are the correlation functions and conditional pair 
distributions that were defined in Section [H] Figure IT51 
shows representative correlation functions between two 
charged hadrons, while Figs. 1131 and 1141 show represen- 
tative conditional yield distributions triggered on neu- 
tral pions (ir°) and charged pions respectively. All three 
correlation sets (Fig. ^] to I14fl show relatively narrow 
peaks centered at A0 = and tt radians. The widths 
of these structures decrease with larger px, which is con- 
sistent with narrowing of the jet cone for increasing p^. 
The fractional area under the jet peak relative to the 
flat underlying background also increases significantly as 
function of associated particle px, indicating increasing 
(di)jet contributions to the correlation function. In par- 
ticular, Fig. El shows that for events where there is a 
high px trigger, a large fraction of the low px (as low as 
0.4—1 GeV/c) particles are coming from the dijet frag- 
mentation, and the jet contribution dominates at px > 2 
GeV/c. Events tagged with a high px jet are much harder 
than a typical minimum bias event. 

We characterize the jet correlations shown in Figs. 1121 
to 1141 by assuming that there are only two contributions 
to the correlation function - (di)jet correlations and an 
isotropic underlying event. This scenario can then be 
expressed as: 
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TABLE III: Summary of the systematic errors on the widths and jr, (sin 2 (4>jj)} ■ 



Error source 



tracking cuts, pair cuts < 4% 

assumptions used in formula < 5% 

S/B correction (n° only) 1 - 3% 



TABLE IV: Summary of the systematic errors on the Conditional Yields for — h analysis. 

Error source < 4 GeV/c 4 — 6 GeV/c 

quality cuts < 4% < 4% 

background correction 3% 30% 

error on single particle yields 10% 10% 



TABLE V: Summary of the systematic errors on the Conditional Yields for 7r° — h analysis 
Pair cuts 3% 

Single Particle 

normalization 0.070 
£smglc p smearing (reso+scale) 3% 



near-side yield 1% 



S/B PT,assoc (GeV/c) <2 2-3 >3 

far-side yield 5% 2% 1% 



TABLE VI: Summary of the systematic errors on the Conditional Yields for ty — h analysis. 



Single Particle 

^single 



normalization 
p smearing (reso+scale) 
trigger pion background 
centrality dependent part 



6.5% 

3% 
5% 
5% 



Yield Extraction 



PT.assoc (GeV/c) 



< 1 



1 - 2 



2-3 



pair cuts 
near-side yield 
far-side yield 
error on the fit 



1% 
20% 

10 - 20% 



1% 
10% 

6% 



2% 
6% 



6% 



4% 



3-4 



3% 
6% 

4% 



4% 
6% 

4% 



TABLE VII: Summary of the systematic errors on the p ou t distribution for 7T — h analysis. 



Pout (GeV/c) 


< 0.5 


0.5- 1 


1 - 2 


2 - 2.5 


yield extraction (near) 


8% 


15% 


20% 


20% 


yield extraction (far) 


8% 


15% 


20% 


30% 


other errors 


10.6% 


10.6% 


10.6% 


10.6% 



C(A0) = A (l + J(A</>)) 



(36) 



where A D denotes the isotropic background and J(A<j>) 
is the jet-function. Approximating the jet-function as the 
sum of two Gaussians, we fit the correlations with; 



Here, \n,f are the normalized Gaussian areas and <jn,f 
are the Gaussian widths for the near and far-side jets re- 
spectively. For the pair distribution functions we fit with 
the same shaped function, but with a different normal- 
ization (Eq. I26fl as outlined in section ITlI Bl 

Figure ^| shows the associated-pr dependence of the 
extracted widths for both the near- and far-side peaks 2 



C(Ao» = A (l 



A 



JV 



27TCTJV 



A, 



2tt(tf 



)• 

(37) 



The results here are not sensitive to the slightly different range 
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25 < PT.trig < 4 GeV/C 4 < P T ,.rig 6 GeV/C 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

A $ (rad) A <|> (rad) 



FIG. 12: Assorted ft — ft correlation functions from d + 
Au collisions for centrality 0-80% and several pr cuts. The 
trigger pr range is 2.5 < pr.trfg < 4.0 GeV/c (a - b) and 4.0 < 
PT,tri g < 6.0 (c - d) respectively, while the associated hadron 
is in the range 1.5 < p T < 2.0 GeV/c or 2.0 < p T < 2.5 
GeV/c. The correlations are for the centrality class — 80%. 
The dashed line represents a fit to the correlation function 
using Eq. E3 



from the charged-hadron correlation functions with the 
trigger range for the charged hadron being 3 — 5 GeV/c. 
The data are tabulated in Table IVlTll 



TABLE VIII: Near and far-side widths as a function of 
f>T,assoc for charged hadron triggers (3 — 5 GeV/c) and as- 
sociated charged hadrons from d + Au collisions. 



(jpT, assoc ) 


&near 


<J far 


(GeV/c) 


(rad) 


(rad) 


0.59 


0.411 ±0.055 


0.89 ±0.28 


0.83 


0.395 ± 0.039 


0.807 ±0.128 


1.2 


0.364 ± 0.032 


0.636 ± 0.079 


1.7 


0.291 ±0.023 


0.688 ±0.103 


2.2 


0.246 ±0.019 


0.637 ±0.146 


2.7 


0.236 ± 0.023 


0.415 ±0.114 



In Fig. 1161 we present the same quantities from the 



-©- 
< 
-a 



3 



1 

0.5 

0.2 
0.1 



7i* -h* -d + Au *p + p 

0.4 < p Tfl9soc < 1.0 GeV/c 
,'\ ° Am 



2.0<p Ta5Soc <3.0 GeV/c 



0.6 
0.4 
0.2 

0.2 
0.1 



2 4 2 4 



1.0<p Tflssnc <2.0 GeV/c 



•-o o-— V ^ % "> 



4* 



3.0 < p Taasoc < 5.0 GeV/c 



A(|)(rad) 



A(|)(rad) 



FIG. 13: Fully corrected assorted charged pion-hadron con- 
ditional pair distributions for d + Au collisions centrality 
— 88% and p ± p collisions. The trigger ir s are within 
5 < PT,tri g < 10 GeV/c and are correlated with hadrons with 
PT.assoJ 0.4 - 1.0 GeV/c, 1.0 - 2.0 GeV/c, 2.0 - 3.0 GeV/c 
and 3.0 — 5.0 GeV/c (from top to bottom and left to right). 
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FIG. 14: Fully corrected assorted 7r°-hadron conditional pair 
distributions for d ± Au collisions centrality — 88% and 
p + p collisions. The trigger 7r°'s are within 5 < pT.trig < 
10 GeV/c and are correlated with hadrons with pT,assoc of 
(upper-left) 1.5 — 2 GeV/c, (upper-right) 2 — 3 GeV/c, (lower 
left) 3-4 GeV/c, and (lower right) 4-5 GeV/c. 



high-pT identified pion correlations, where there is ex- 
cellent agreement between the 7T° and charged-pion data 
sets. For both types of identified pions the trigger px 
range is 5 — 10 GeV/c and these data are tabulated in 
Tables llXl and H 

The far-side widths shown in Figs . 1 1 51 and 1 1 61 are larger 
than the near-side widths, as expected, since the far-side 
structure is a convolution of two jet fragmentations as 
well as any kx of the scattered partons. The widths of 
the correlation functions also steadily decrease as a func- 
tion of pT.assoc as expected from (di)jet fragmentation. 
For completeness we also tabulate the near- and far-side 
widths extracted as a function of PT,trig for identified pi- 
ons. These data are tabulated in Tables 1x71 and IXIII 

Although PHENIX single particle acceptance is lim- 
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Centrality (%) = 0-80 
(a) near side (filled) 
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FIG. 15: Near (a) and far-side (b) widths as a function of 
PT.assoc for charged hadron correlations from minimum-bias d 
± Au collisions (see text). Bars are statistical errors. 



TABLE IX: Near and far-side widths as a function of pT.assoc 
for charged pion triggers (5 — 10 GeV/c) and associated 
charged hadrons from minimum-bias d ± Au collisions. 



(pT.assoc) 




& far 


(GeV/c) 


(rad) 


(rad) 


0.50 


0.440 ± 0.044 


0.651 ±0.052 


0.70 


0.391 ±0.026 


0.587 ±0.039 


0.90 


0.331 ±0.023 


0.613 ± 0.044 


1.23 


0.271 ±0.010 


0.517 ±0.024 


1.75 


0.210 ±0.008 


0.433 ± 0.022 


2.24 


0.193 ±0.009 


0.372 ± 0.023 


2.73 


0.165 ±0.007 


0.317 ±0.020 


3.44 


0.135 ±0.006 


0.307 ±0.020 


4.42 


0.128 ±0.008 


0.287 ±0.023 



ited to \rf\ < 0.35, it can sample jet pairs in twice as 
large of a window in Ar] (\Arj\ < 0.7) with varying pair 
efficiency. Similar to azimuthal correlation, the pair effi- 
ciency in Arj can be estimated via mixed-events and can 
subsequently be divided out (Eq. I21|) . Assuming that 
the underlying event is flat 3 in |A?7| < 0.7, we fix the 
background level to be equal to that in azimuthal corre- 
lation function Eq. 1361 and subsequently extract the jet 
distribution as a function of Ar). In Fig. 117b . we com- 
pare the near-side jet shape in Acf> and Arj in the angular 
range of | A<fi, Arj\ < 0.7 for the — correlation with 
1.0 < PT,assoc < 2.0 GeV/c. There is no significant dif- 
ference in jet shape between A77 and A<p and the widths 



3 rj depe nden ce of the single particle yield is very weak in < 
|r?| < 1 l4.'l . Thus the underlying pair distribution in |A?j| < 0.7 
is almost flat. 



TABLE X: Near and far-side widths as a function of pT.assoc 
for neutral pion triggers (5 — 10 GeV/c) and associated 
charged hadrons from minimum-bias d ± Au collisions. 



(PT.assoc) 




& far 


(GeV/c) 


(rad) 


(rad) 


1.21 


0.284 ±0.011 


0.494 ± 0.022 


1.71 


0.227 ±0.007 


0.410 ±0.019 


2.37 


0.193 ±0.005 


0.380 ±0.015 


3.39 


0.177 ±0.006 


0.322 ± 0.020 


4.41 


0.130 ±0.007 


0.315 ±0.026 
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FIG. 16: (Color online) a) near-side width, b) far-side width 
as a function of pr.assoc for a charged pion (closed symbols) 
and neutral pion (open symbols) triggers from the 
range of 5 — 10 GeV/c in minimum- bias d ± Au collisions 
(see text). Bars are statistical errors. 



are consistent in both directions. Figure IT7b shows the 
far-side jet shape in Ar), the associated pair distribution 
is flat within ±10%. 

Figure 118b shows the comparison of the near-side jet 
width in A<j> and A77 from d ± Au. There is overall very 
good agreement between the two data sets. However, 
the width in A77 is systematically lower than that in A(j> 
at small pr.assoc- This is due to the fact that the the 
underlying background is not completely flat in A77, but 
varies by up to 10% in |A?j| < 0.7. Thus the procedure of 
dividing by the mixed-event distribution (Eq. |2J intro- 
duces some distortion of the jet shape at large Ar], and 
consequently leads to a slightly different value for the jet 
width. In fact for p + p collisions, Figure HHb indicates a 
similar discrepancy between A<f> and A77 at small PT,assoc 
for p + p collisions. Thus this deviation is not likely due 
to the cold medium effect in d ± Au. 

We extract not only the widths of the jet-structures 
but also the conditional yields of how many hadrons are 
in the near-side and far-side structures for each high- 
Pt trigger. The conditional yield defined in Eq. 1141 can 
be obtained from either a correlation function or condi- 
tional pair distribution, both of which produce identical 
results. For the conditional pair distributions the condi- 
tional yield is directly extracted from the fit parameters 
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TABLE XI: Near and far-side widths as a function of PT,tri g 
for charged pion triggers and associated charged hadrons (2 — 
4.5 GeV/c) from minimum-bias d + Au collisions. 



(PT.trig) 

(GeV/c) 


&near 

(rad) 


(rad) 


5.44 


0.176 ±0.008 


0.393 ± 0.030 


6.31 


0.165 ±0.007 


0.342 ± 0.020 


7.27 


0.162 ±0.007 


0.322 ± 0.022 


8.60 


0.157 ±0.008 


0.301 ±0.019 


10.6 


0.149 ±0.020 


0.231 ±0.039 


13.2 


0.177 ±0.019 


0.329 ± 0.042 


TABLE XII: Near and far-side widths as a function of PT.trfg 
for neutral pion triggers and associated charged hadrons (2.5— 
5 GeV/c) from minimum-bias d ± Au collisions. 


(PT.trig) 

(GeV/c) 


(rad) 


Pfar 

(rad) 


5.39 


0.207 ±0.007 


0.409 ± 0.025 


6.40 


0.151 ±0.008 


0.304 ± 0.034 


7.66 


0.144 ±0.012 


0.295 ± 0.035 



(Eg. I26f) . while for correlation functions several normal- 
ization factors need to be applied to obtain the per trigger 
yield [44|,|45(, as described below. 

For correlation functions, it is convenient to define 
the fraction of jet-correlated particle pairs per event, 
" Jctp °' r . Following the basic ansatz outlined in Ea. 1361 
the fraction of jet-correlated particle pairs is obtained by 
summing the jet function over all bins in A(f> and dividing 
by the total sum of the correlation function 




dN/A4>, a = 0.250 ±0.007 
dN/Ar|, a = 0.234 ±0.007 
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FIG. 17: The jet shapes in A(f> and Arj from ir^ — h corre- 
lation with 5 < PT.trig < 10 GeV/c and 1.0 < pT,a SS oc < 2.0 
GeV/c from minimum- bias d ± Au collisions; a) The near- 
side jet shape in A(j> (solid circles) and Arj (open boxes), b) 
The far-side jet shape in Ar/. 



"'jet pair 
^total pair 



EA,J(A0) 
£C(A0) ' 



(38) 



Such pair fractions are shown as a function of pr.assoc 
for a trigger hadron of 3.0 < pr < 5.0 GeV/c and a 
centrality selection of — 80% in Fig. ED The results, 
shown for both the near and far-side jets, indicate an 
increase in the average fraction of jet-correlated particle 
pairs with px as might be expected if jet fragmentation 
becomes the dominant particle production mechanism as 
Pt is increased. 

The pair-fraction is multiplied by the ratio — 



fc jct pair 



fc jet pair 



^trig^assoc ^total pair ^trig^assoc 



(39) 



where n pa i rs denotes the average number of detected par- 
ticle pairs per event and n tr ig, and n ass0 c are the detected 
single particle yields per event for trigger and associated 
particles respectively. This gives the average number 
of jet-correlated pairs per event over the combinatoric 



background 



The conditional per-trigger yield, 



' Jctpa ' r , is obtained via multiplication by the efficiency 

3soc° rr ') f° r t ne selected 



corrected single particle yield (n|f £ 
associated pt bin of interest; 



"jet 1 



'-jctpair 



^eff — corr. 



Tlt r igTl a 



(40) 



The per trigger yields for hadron triggers (found using 
Eq. I40fl are corrected for the azimuthal acceptance and 
tracking efficiency but are reported within the PHENIX rj 
acceptance for the central arms, i.e. no R(Ar)) corretions 
is applied to the hadron-triggered conditional yields. 

Figure 1201 plots the near- and far-side invariant con- 
ditional yields extracted via Eq. 0U| for different trigger 
Pt selections as indicated. An approximate exponential 
decrease with px is observed, i.e. there are more low- 
Pt particles associated with each high-pT trigger hadron. 
The data are tabulated in Tables IxTTTI and IxTVl 

In Fig. [2] the conditional yields for identified pion 
triggers are plotted as function of pr.assoc for both near- 
side correlation and far-side correlation. For this high-py 
data, the conditional yields are extracted from the fits to 
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FIG. 18: The comparison of jet width as function of PT.assoc 
in A<j) (solid circles) and Ar/ (open boxes) from ir — h cor- 
relation with 5 < PT.trig < 10 GeV/c. a) results for d ± Au. 
b) results for p + p. Bars are statistical errors. 
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FIG. 19: Average jet pair fraction per event as a function of 
PT.assoc- Results are shown for the trigger hadron selection 
3.0 < pt < 5.0 GeV/c from d + Au collisions and a centrality 
of — 80%. Bars are statistical errors. 



the data in Figs . 1 1 31 and 1 1 41 using Eq. [2H1 then corrected 
for pair efficiency in A77 using Eq. [21 and Eq. [531 The 
conditional yields are tabulated in Tables IXVl and IXVll 

The agreement between the two pion-triggered data 
sets is good, which indicates that the jet fragmentation 
function is independent of whether a neutral pion or a 
charged pion trigger is used. The difference in the mag- 
nitudes of far-side and near-side yield reflect the fact 
that the far-side correlations measures a hadron trig- 
gered effective fragmentation while the near-side corre- 
lation measures dihadron fragmentation. 

The conditional yields presented in FigOAs. [201 an d 
l^can be considered as the basic information, while the 
near and far dN/dxE distributions in Section TlV Bl have 
a closer relationship to parton fragmentation functions 
as was described in Section ITT1 

Since multiple scattering should increase with central- 



TABLE XIII: Near and far-side conditional yields as a func- 
tion of p-r.assoc for charged hadron triggers (2.5 — 4 GeV/c) 
and associated charged hadrons from d + Au collisions. 



(PT.assoc) 

(GeV/c) 


dN / dpTnear 


dN/dp T far 


0.592 


0.327 ±0.092 


0.383 ±0.182 


0.831 


0.307 ± 0.045 


0.339 ± 0.079 


1.190 


0.174 ±0.018 


0.158 ±0.024 


1.702 


0.081 ± 0.009 


0.066 ±0.014 


2.205 


0.042 ± 0.006 


0.028 ± 0.009 


TABLE XIV: Near and far-side conditional yields as a func- 
tion of pr.assoc for charged hadron triggers (4 — 6 GeV/c) and 
associated charged hadrons from d + Au collisions. 


(pT,assoc) 

(GeV/c) 


dN / dprnear 


dN/dprfar 


0.831 


4.437 ± 1.040 


6.031 ±2.010 


1.200 


2.725 ±0.506 


2.051 ±0.562 


1.700 


1.907 ±0.278 


2.046 ± 0.447 


2.210 


0.819 ±0.152 


0.804 ± 0.244 


2.931 


0.497 ±0.107 


0.258 ±0.061 



ity, we examine whether these jet characteristics exhibit 
any centrality dependence. Figure |23 reports on widths 
and conditional yields in the PHENIX r\ acceptance for 
dihadron correlations with a trigger hadron in the range 
2.5 < PT,trig < 4 GeV/c. The data are tabulated in Ta- 
ble lXVIli and lX Villi Centrality dependent widths for the 
7T°— h correlations are shown in Fig. |22|for a mean 7r° pt 
of approximately 5.4 GeV/c. The data for the widths are 
tabulated in Table 1X1X1 while the yields can be found in 
Table IXXl 

Neither the widths nor the per-trigger yields change 
significantly with centrality, indicating that the influence 
of multiple scattering on jet properties is small in this 
region. This work is extended in Section llV Dl where we 
present the centrality dependence of various jet-structure 
observables. 



B. Jet Properties in Minimum Bias d ± Au 
Collisions 



From the angular widths and yields in the previous 
section, we calculate the following quantities that char- 
acterize the jet structures: transverse momentum of 
hadrons with respect to the jet (jr)> the dijet acopla- 
narity (sm 2 ((f>jj)) , and the dN/dxE distributions. These 
quantities are first presented for minimum-bias d ± 
Au collisions, which have the highest statistical preci- 
sion, and are then compared with results from p ± p in 
Section ITVCl 

Figure |2] shows the compilation of \J values ex- 
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FIG. 20: (Color online) Per trigger yield (a) and invariant 
conditional yield (b) as a function of pr.assoc for trigger hadron 
ranges of 2.5 < p T < 4.0 GeV/c and 4.0 < p T < 6.0 GeV/c, 
respectively. The closed points are the near-side yields and 
the open points are the far-side yields. The centrality range 
is — 80% in d + Au collisions. The yields are corrected for 
efficiency and are reported in the PHENIX eta acceptance. 
Bars are statistical errors. 



tracted from n° — h correlations (open crosses), ^ — h 
correlations (filled circles) , and assorted pt and fixed pt 
correlation results for h— h at low px shown in solid boxes 
and solid stars, respectively. The yj (j T ) values were cal- 
culated using Eq. The systematic errors are mainly 
due to the uncertainties of Eq. [7| and the fitting proce- 
dure, which are about 5% (independent of pt) and are 
approximately the same for all four analysis. The U (j' T ) 
values for — h and tt° — h indicate a steady increase 
at p T below 2 GeV/c followed by a saturation around 
560 - 640 MeV/c at p T > 2 GeVJc. The initial increase 
is due to the 'seagull effect' 0,0. The approximately 
constant behavior of y/Jjx) above 2 GeV/c is consistent 
with the scaling behavior of the fragmentation functions. 
The \J (j' T ) results from the h — h correlation analyses 
have a similar increase and saturation behavior. They 
seem to reach a slightly higher plateau value at a lower 
Pt (around 1.5 GeV/c). 



TABLE XV: Near and far-side conditional yields as a func- 
tion of p-r.assoc for charged pion triggers (5 — 10 GeV/c) and 
associated charged hadrons from minimum-bias d + Au col- 
lisions. 



(pT.assoc) 

(GeV/c) 


dN/dprnear 


dN/dprfar 


0.5 


1.57 ±0.083 


2.54 ±0.15 


0.7 


0.911 ± 0.049 


1.53 ±0.13 


0.9 


0.574 ± 0.031 


1.00 ±0.087 


1.1 


0.542 ± 0.032 


0.727 ±0.068 


1.3 


0.456 ± 0.026 


0.801 ± 0.058 


1.5 


0.351 ±0.022 


0.451 ± 0.044 


1.7 


0.303 ±0.018 


0.365 ± 0.038 


1.9 


0.235 ±0.015 


0.327 ±0.031 


2.1 


0.172 ±0.012 


0.222 ± 0.025 


2.3 


0.135 ±0.010 


0.203 ± 0.022 


2.5 


0.108 ±0.008 


0.162 ± 0.018 


2.7 


0.0905 ± 0.0075 


0.145 ±0.017 


2.9 


0.0742 ± 0.0064 


0.107 ±0.013 


3.1 


0.0645 ± 0.0059 


0.070 ±0.011 


3.3 


0.0490 ± 0.0052 


0.0819 ±0.0113 


3.5 


0.0473 ± 0.0047 


0.0647 ± 0.0097 


3.7 


0.0439 ± 0.0045 


0.0636 ± 0.0084 


3.9 


0.0367 ± 0.0040 


0.0495 ± 0.0075 


4.1 


0.0281 ± 0.0034 


0.0327 ± 0.0064 


4.3 


0.0297 ± 0.0034 


0.0446 ± 0.0068 


4.5 


0.0256 ± 0.0031 


0.0238 ± 0.0048 


4.7 


0.0192 ±0.0027 


0.0397 ± 0.0061 


4.9 


0.0112 ±0.0021 


0.0137 ±0.0036 


TABLE XVI: Near and far-side conditional yields as a func- 
tion of PT,assoc for neutral pion triggers (5 — 10 GeV/c) and 
associated charged hadrons from minimum-bias d + Au col- 
lisions. 


(pT,assoc) 

(GeV/c) 


dN/dpTnear 


dN/dprfar 


1.21 


0.545 ± 0.035 


0.718 ±0.076 


1.71 


0.289 ±0.014 


0.389 ± 0.030 


2.21 


0.236 ± 0.009 


0.203 ±0.018 


2.72 


0.102 ±0.006 


0.155 ± 0.012 


3.22 


0.0724 ± 0.0049 


0.090 ± 0.012 


3.73 


0.0448 ± 0.0039 


0.0560 ± 0.0071 


4.23 


0.0308 ± 0.0029 


0.0481 ± 0.0058 


4.72 


0.0152 ±0.0017 


0.0495 ± 0.0060 



A combined fit based on U (J'i,) data points at pr > 
2 GeV/c gives a plateau value of \J (j T ) — 0.64 ± 
0.02(stat) ± 0.04(sys) GeV/c for minimum bias d + Au 
collisions. 

A key quantity that provides information on multi- 
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FIG. 21: Conditional yield as function of pT.assoc for near- 
side (left panel) and far-side (right panel) correlation from 
minimum-bias d + Au collisions. Bars are statistical errors. 
The shaded boxes represent the total systematic errors on 
each point. 



TABLE XVII: Near and far-side widths as a function of N co ;z 
for charged hadron triggers (2.5 — 4 GeV/c) and associated 
charged hadrons (1-2.5 GeV/c) from d ± Au collisions. 



Centrality 


<Jncar{rad) 


af ar (rad) 


0-20 % 


0.351 ±0.021 


0.669 ± 0.072 


20-40 % 


0.364 ± 0.024 


0.670 ± 0.068 


40-80 % 


0.325 ±0.018 


0.628 ± 0.070 



pie scattering in the cold nuclear medium is (sin 2 (<pjj)), 
where <pjj is the azimuthal angle between the jet axes. 
As described in Section [H] we calculate (sin 2 ^-) from 
the experimental values of the near- and far-side widths 
(Eq. 1 1 3|> - Figure [23 shows (sin 2 0jj) as function of 
PT.assoc for high-pT pion triggers. We observe that the 
RMS of the sine of the angle between the jet axes, 
(sin 2 <j>jj), decreases as higher-p^ associated particles are 
selected. There is good agreement between the two data 
sets. Figure plots (sin 2 <j>jA as a function of px, trig, 
where a similar decrease with px is observed. 

In the next section we will calculate the quadrature 
difference (sin 2 <j)jj) between d±Au and p + p collisions, 
and use that to quantify the affect of additional final- 
state scattering in <i±Au collisions. 

In Section[n]we defined the near-side and far-side p ou t- 
With this observable, it is possible to move beyond cal- 
culating means or RMS values and hence the p out dis- 
tribution potentially carries more information about the 
dijet acoplanarity. The measured p ou t distributions for 
h ± are shown in Fig. [23 for the near-side and far- 
side. The far-side p on t bas a broader distribution than 
the near-side p ou t, reflecting the fact that kx is larger 
than jx- The p ou t distributions have a power law tail, 
possibly due to strong radiative processes driving large 
values of p ou t ■ 

In Section IIV Al we reported the yields of associated 
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FIG. 22: Near- and far-side widths as well as conditional 
yields in the PHENIX pseudorapidity acceptance for several 
centrality selections from d + Au collisions. Results are shown 
for the trigger hadron selection 2.5 < pr < 4.0 GeV/c and 
the associated pr range of 1 < pr < 2.5 GeV/c. Bars are 
statistical errors. 



TABLE XVIII: Near and far-side conditional yields as a func- 
tion of N co n for charged hadron triggers (2.5 — 4 GeV/c) and 
associated charged hadrons (1-2.5 GeV/c) from d + Au colli- 
sions. 



Centrality 


dN/ dpxnear 


dN/dpxfar 


0-20 % 


0.150 ±0.021 


0.113 ±0.030 


20-40 % 


0.182 ± 0.023 


0.118 ±0.033 


40-80 % 


0.151 ± 0.014 


0.102 ±0.023 



hadrons per trigger particle, or the conditional yield. A 
more comprehensive way of quantifying the fragmenta- 
tion function is to plot the conditional yields as a func- 
tion of xe- This is shown in Fig. [2H1 for ^ — h (solid 
circles) and 7r° — h (open crosses). 

Previously, in ISR experiments, the slope of the 
xe distribution has been determined to be around 5.3 
(GeV/c)- 1 HJ. In Fig. E3 the conditional yields as 
a function of i:^ are plotted for trigger px range of 
5 — 6 GeV/c. In order to compare data with the pre- 
vious ISR results, we have determined the exponential 
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FIG. 23: Near- and far-side widths for several centrality se- 
lections from d ± Au collisions. Results are shown for 7T° — h 
correlations with the trigger 7r° 5 < pr < 10 GeV/c and the 
associated pr range of 2 < pr < 3 GeV/c. Bars are statistical 
errors. 



TABLE XIX: Near and far-side widths as a function of cen- 
trality for neutral pion triggers (5 — 10 GeV/c) and associated 
charged hadrons(2-3 GeV/c) from d + Au collisions. 



Centrality 


a„ear(rad) 


a far (rad) 


0-20 % 


0.199 ±0.009 


0.387 ± 0.024 


20-40 % 


0.195 ±0.009 


0.401 ±0.031 


40-88 % 


0.190 ±0.008 


0.376 ± 0.024 



TABLE XX: Near and far-side conditional yields as a function 
of centrality from d ± Au collisions for neutral pion triggers 
(5 — 10 GeV/c) and associated charged hadrons(2-3 GeV/c). 



Centrality 


dN/dprnear 


dN/dpTfar 


0-20 % 


0.0816 ±0.0037 


0.116 ±0.008 


20-40 % 


0.0947 ± 0.0045 


0.141 ± 0.011 


40-88 % 


0.0967 ± 0.0043 


0.144 ±0.009 
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FIG. 24: (Color online) the extracted yj <j|,) for minimum 
bias d ± Au collisions from all four dihadron correlations. The 
trigger pr ranges are 3 — 5 GeV/c, 5—10 GeV/c and 5 — 10 
GeV/c for h ± — ft* (filled boxes), n° — /i* (open crosses) and 

— fe (filled circles) assorted-pr correlations. The h — h 
fixed-pr correlation is shown by filled stars. The statistical 
and systematic errors are combined for the assorted-pr cor- 
relations, while the fixed-pT correlation shows the systematic 
error separately as shaded boxes. 



inverse slope 0.3 < xe < 0.7 and obtain the inverse slope 
parameter of 6.0 ± 0.3 (GeV/c) -1 in the near-side and 
7.1 ±0.5 (GeV/c) -1 in the far-side. The near-side Xe in- 
verse slope is smaller than that for the far-side, reflecting 
the difference between dihadron fragmentation and sin- 
gle hadron fragmentation. By requiring a trigger particle 
on the near-side, one reduces the jet energy available for 
fragmenting to the second hadrons and consequently a 
smaller inverse slope for the near-side. Note, however, 
that the xe distributions do not have pure exponential 
shape, and the fitted inverse slope is sensitive to the fit- 
ting ranges. 

It is well known that fragmentation functions D(z) ap- 
proximately scale in e + e~ or p + p collisions, i.e. are 
independent of jet energy. To check whether this is still 
true in d ± Au collisions, we plot in Fig. [30 the condi- 
tional yields as a function of Xe for different ranges of 
trigger px from n^ 1 — ft, ± correlations. All curves fall on 
top of each other, indicating a universal behavior of the 
jet fragmentation function. 

At lower px we have xe distributions from the h- 



h correlations. In Fig. |^ far-side conditional yields 
as obtained from charged hadron correlation functions 
are plotted versus (xe). Here, (xe) has been calcu- 
lated from the (pr.trig), (pr.assoc) and extracted angu- 
lar widths. Results are shown for two trigger px ranges, 
2.5 < pr.trig < 4 GeV/c and 4 < ptmi S < 6 GeV/c, 
respectively. The dashed lines represent exponential fits 
to the data. The slopes extracted are 6.3 ± 1.2 for the 
lower trigger px range and 6.1 ± 0.8 for the higher trig- 
ger px window, respectively. Within the statistics of the 
charged hadron dataset we do not observe a strong sen- 
sitivity of the slope of the xe distributions to the trigger 
p T . 

A direct way of quantifying the scale dependence of the 
xe distribution is to plot the far-side conditional yields 
versus pr.trig for a fixed range of xe- This is shown in 
Fig. where the conditional yields are found to be in- 
dependent of pt, trig, i-e. there is no significant deviation 
from scaling. We will quantify any scaling violation in 
this data when we compare these xe distributions from d 
± Au to distributions from p+p collisions in Section FlV CI 
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FIG. 25: (Color online) (sin 2 for minimum bias d + Au 
collisions as function of associated particle pr for ir^ — h 
(filled circles) and 7r° — h (open crosses) correlations where 
the trigger particles have a pr between 5—10 GeV/c. The 
data points are plotted at the mean pT.assoc- Bars are statis- 
tical errors. The boxes represent the total systematic errors 
on each point. 
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FIG. 26: (Color online) (sin 2 <j)jj\ for minimum bias d + Au 
collisions as function of trigger particle pr for tt ± — h (filled 
circles) and ir° — h (open crosses) correlations. The data 
points are plotted at the mean pr.trig- Bars are statistical 
errors. The boxes represent the total systematic errors on 
each point. 



C. Comparison between d + Au and p + p 

As discussed in Section [I] multiple scattering in the 
cold nuclear-medium may broaden the far-side correla- 
tion and possibly modify the fragmentation properties. 
In the previous section we presented the measured jet 
structures from minimum-bias d + Au collisions. In this 
section we compare that data to results from p + p col- 
lisions. The goal is to establish the extent to which the 
nuclear-medium modifies the properties of jets. 

Figure|2Slshows the comparison of the extracted y/TJj) 
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FIG. 27: (Color online) Near-side and far-side p out distri- 
butions for minimum bias d + Au collisions obtained from 
— h correlation. The trigger is 5 < PT.trig < 10 GeV/c, 
the associated particle is 0.5 < pT.assoc < 5.0 GeV/c. Bars 
are statistical errors. The boxes represent the total systematic 
errors on each point. 



10 



HI 

X 
73 



10- 1 



6 <P T , t r ig < 10GeV/C 



•» 



7t°- If 
7ft - hf 



Near side 



V 



+ 



Far side 



7" 



0.2 0.4 0.6 0.E 0.2 0.4 0.6 0.8 



E,near 



E,far 



FIG. 28: Conditional yield as a function of xe for near- 
side (left panel) and far-side (right panel) for tt^ — h (filled 
circles) and ty° — (open crosses) from minimum-bias d + 
Au collisions. Bars are statistical errors. The boxes represent 
the total systematic errors on each point. 



values as function of PT,assoc from d + Au and p + p col- 
lisions. The yj (jj,) values show no change from p + p to 
d + Au collisions. 

Similarly, the (sin 2 {(f>jj)) values shown in Fig. [M] for d 
+ Au are comparable to those from p + p within errors, 
although the values from d + Au collisions are system- 
atically higher for — h . Since there is no strong 
difference between the cZ+Au and p + p results, there is 
little indication for increased multiple-scattering in the 
d+Au final state. 

Any additional radiation can be quantified by calculat- 
ing the point-by-point quadrature difference in (sin 2 <j)jj) 
between d + Au and p + p collisions. This is shown 
in Fig. I35I and this difference is consistent with zero. 
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FIG. 29: (Color online) Conditional yield as a function of xe 
for near-side (left panel) and far-side (right panel) correlations 
for ir —h correlations from minimum-bias d + Au collisions. 
The trigger pions are 5 < f>T,trig < 6 GeV/c, and the black 
lines are fits to an exponential function for 0.3 < xe < 0.7. 
Bars are statistical errors. 
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FIG. 30: (Color online) Conditional yield as a function of xe 
for near-side (left panel) and far-side (right panel) correlation 
for 7T — h correlation for several different trigger pr's for 
minimum-bias d + Au collisions. Bars are statistical errors. 



The average value for ir° — h is A (sin 2 (f>jj) = —0.005 ± 
0.012(stat)±0.003(sys); for ^-h A (sin 2 ^) = 0.011 + 
0.011(stat)±0.010(sys). Combining the two data sets, we 
find A (sin 2 = +0.004 ± 0.008(stat)±0.009(sys). 

Figure shows the comparison of the p ou t distribu- 
tion between central d + Au and p + p from ■n ± — h 
correlation, no apparent differences are observed in both 
the near and far-side. This is consistent with the ob- 
servations that both yj (jj,) and (sin 2 (4>jj)) are similar 
between d + Au and p + p. 

A second set of comparisons between d + Au and p + p 
collisions is the number of hadrons in the near- and far- 
angle jet structures associated with a high-p^ trigger. 
Figure I37I shows the comparison of the conditional yield 
as function of xe and no apparent difference between d 
+ Au and p + p collisions is observed for either the near- 
or far-side. 

In the previous section we examined the level of scaling 
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FIG. 31: (Color online) Conditional yield function 
of (xe) for the far-side h — h correlations. Data is shown 
for two trigger pr ranges, 2.5 < PT,tri g < 4 GeV/c (open 
circles) and 4 < pr,tri g < 6 GeV/c (open boxes), respectively. 
The dashed lines are exponential hts to the data. Bars are 
statistical errors. 
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FIG. 32: Far-side conditional yield as a function of pr.trig for 
different ranges of xe for — h correlation for minimum- 
bias d + Au collisions. Bars are statistical errors. 



violations in the far-side dN / dx e distribution for d + Au 
collisions by plotting different x e ranges as a function of 
PT,trig (Fig. The comparable plot for p + p collisions 
is shown in Fig. EHl For both d + Au and p + p collisions 
the amount of scaling violations, i.e. the dependence of 
dN/dxE on pr.trig, can be quantified by fitting the data 
in each xe range with a straight line as a function of 

PT,trig- 



dN 

dXE 



dXEo 



■Ppt) 



(41) 



The fitted slopes (/3) represent the fractional change in 
dN/dxE per GeV/c and are shown in Fig. For the 
d + Au data (3 is consistent with zero, i.e. there is no 
significant scaling violation across the whole xe range, 
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FIG. 33: (Color online) The comparison of the y(ir) values 
between d + Au (filled circles) and p + p (open circles) for 
n ± — h correlations (top panel) and 7T° — h correlations 
(lower panel). The trigger pion range is 5-10 GeV/c. Bars are 
statistical errors. The systematic errors are given in Table liTTl 



while there may be a slight scaling violation at high xe 
for the p + p data. On a point-by-point basis there is no 
systematic difference between the d + Au and p+p data. 

Taken as a whole, all the results presented in this sec- 
tion indicate that there is no significant change in jet 
fragmentation between d + Au and p + p collisions due 
to the presence of the cold nuclear- medium. In addition 
there is no strong evidence for an increase in (sin 2 (</>jj)) 
due to multiple scattering in the Au nucleus. Using the 
minimum-bias d + Au data has the advantage of the 
highest statistical precision. In the next section we ex- 
amine whether there is any change in jet-structures as a 
function of collision centrality in d + Au, i.e. we split 
the statistics into a few centrality classes to increase the 
lever arm of the nuclear-thickness function. 



D. Centrality Dependence 

As discussed in Section (sin 2 ^^)) is expected to 
increase as d + Au collisions become more central due to 
increased multiple scattering. Models of multiple scat- 
tering |12| predict that the increase in (fcy) is propor- 
tional to Ta(6), the nuclear thickness function. We are 
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FIG. 34: (Color online) The comparison of the (sin 2 (<^-)} 
values between d + Au (filled circles) and p + p (open circles) 
for 7T — h correlations (top panel) and it — h correlations 
(lower panel). The associated hadron range is 2 < p-r.assoc < 
4.5 GeV/c for the charged pion triggers and 2.5 < p-r.assoc < 5 
GeV/c for the neutral pions. Bars are statistical errors. The 
boxes represent the total systematic errors on each point. 



not aware of predictions of how (sin 2 ((f>jj)} will change 
with centrality, but note that any increase in (k^) wn ^ 
also increase ^sin (4>jj)')- 

To probe this physics, we have measured angular cor- 
relations in three centrality bins for d + Au collisions 
(0-20%, 20-40%, and 40-88%) to extract angular-widths 
of the jet-structures and hence a/ (j 2 ,) and (sin 2 (<^-j)). 
Figure |4*UI shows the independent data sets of (sin 2 (^)) 
including results from p + p collisions as well as the three 
centrality classes from d + Au collisions. 

All the (sin 2 ^^)) data in Fig.2Dlh ave been simulta- 
neously fit with the following linear equation in Ta(6) 

(sin (4>jj)) = (sin (<fijj ))o(l + a frac T A (b)) (42) 

The slope parameter, af rac , is assumed to be common to 
all data sets, while the pre-factors, (sin 2 '{4>jj))a, depend 
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FIG. 35: (Color online) Quadrature difference between 
minimum-bias d + Au and p + p <^sin 2 <j!>jj^ values. Closed 
circles are 71" — h values and the open circles are ty° — h val- 
ues. Bars are statistical errors. The boxes represent the total 
systematic errors on each point. 
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FIG. 36: (Color online) The comparison of the p ou t distri- 
bution at the near-side (left panel) and far-side (right panel) 
between central d + Au collisions and p + p collisions. Re- 
sults are obtained for 7T — h correlations with the associated 
hadron range 0.5 < pr.assoc < 5 GeV/c and trigger pion range 
of 5 < PT.trig < 10 GeV/ c. Bars are statistical errors. 
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FIG. 37: (Color online) The comparison of the xe distribu- 
tion from — h correlation at the near-side (left panel) and 
far-side (right panel) between minimum bias d + Au collisions 
(filled circles) and p + p collisions (open circles) . The trigger 
are from 5 — 10 GeV/c. Bars are statistical errors. The 
boxes represent the total systematic errors on each point. 



10 



UJ 

X 
T3 



10" 



0.06 < 
0.12 < 
0.18 « 
0.24 < 
0.30 < 
0.36 < 
0.42 < 
0.48 < 



- < 0.12 
H < 0.18 
I < 0.24 

: < 0.30 

I < 0.36 
I < 0.42 
= < 0.48 
I < 0.54 



71* -h± 



10 



PT,trig( GeV/C ) 



FIG. 38: Far-side conditional yield as function of pr.trig for 
different ranges of i£ from p + p collisions, triggers are 7r 
from 5 — 10 GeV/c. Bars are statistical errors. 



Vitev's predicted effect compared to our experimental 
results, we recast Eq. 031into the same form as Ea. 1421 



on the pr of the trigger and associated particles. The ex- 
tracted slope af rac = —0.01 ± 0.40 mb with chi-squared 
per degree-of-freedom, x 2 / 1 ' — 27/22. This slope is con- 
sistent with zero, i.e. we do not observe any significant 
increase in (sin 2 (4>jj)) with centrality. 

This can be compared to predictions from Hwa and 
Wang 16] who assume no increase in fcx with centrality 
to reproduce the Cronin effect data at RHIC and Qiu 
and Vitev 0] who calculate 



(k Tl 



dijet 



2(fc;)i:ac 



0.72 



Tati 



■/bias 



-T A {b) (43) 



0.72 



2{fcj*) vacJ- 1 Arainbias 



T A (b)) (44) 



Hence their prediction for the fractional increase in (fc^) 
with Ta(6) is from 0.51 to 0.72 depending on the range 
0.25 < (M>) vac < 0.35 (GeV/c) 2 suggested by Qiu and 
Vitev [lj] 4 - Though the predicted fractional increase is 
of a different quantity, it should provide an estimate of 
the magnitude of the fractional increase in sin 2 (4>jj))- 
The prediction is slightly larger than one standard devi- 
ation (statistical) from our experimental result, a/ rac = 



To gain some insight into the magnitude of Qiu and 4 For 
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FIG. 39: The fitted fractional change in dN/dxE per unit 
PT,tri g (/3 in Eq. 1411 of the far-side conditional yield for differ- 
ent ranges of xe- The p + p data are shown as open squares, 
while the d + Au data are shown as filled squares. Bars are 
statistical errors. 



—0.01 ± 0.40 mb. If the measured value for < (k^) vac at 
RHIC turns out to be larger than assumed by Qiu and 
Vitev [14| then the predicted fractional increase in Eg. 1441 
will be smaller. 

Barnafoldi et at [15| have also predicted the increase 
in kx due to multiple scattering in d + Au collisions 
at RHIC. They calculate that (fcf.) increases by C — 
0.35(GeV/c) 2 per collision up to the first four collisions, 
then it saturates. Their prediction is 



(k 2 T ) 

(4) 



= (4)0(1 + 
= (4)0(1 + 



0.35 
0.35 



x (40 x T A )) , T A < 0.1 



x (40 x 0.1)) , T A > 0.1 



(45) 



Barnafoldi et at do not provide values for (fc^)o, how- 
ever if we use the range 0.25 < (fef,)o < 0.35 (GeV/c) 2 
suggested by Qiu and Vitev then the fractional in- 
crease with T A is 40 to 56 for T A < 0.1 followed by no 
further increase. This rapid increase is not observable in 
our data set because the model saturates already in the 
most peripheral d + Au bin where (T A ) =0.11 mb -1 . 

As discussed in Section H] inelastic scattering of the 
hard parton in the cold-medium may also increase the 
conditional yields (CY) of hadrons that are associated 
with a high-px trigger. Figure 1411 shows the centrality 
dependence of the extracted CY, together with CY from 
p+p collisions. The difference can be better illustrated by 
taking the ratio (d + Au/p+p) of the per trigger yield as 
shown in Fig. 021 There is a possible increase in near-side 
particle yield for px < 1 GeV/c in the d + Au collisions. 
In other momentum ranges there is no consistent differ- 
ence between yields in d + Au and p+p collisions. The 
recombination model of Ref. (2^ predicted a factor of 
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FIG. 40: The measured values of (sin 2 (<f>jj)} are shown as a 
function of the nuclear thickness function TaQi). The left- 
most points are from p + p collisions and the other points are 
from the three centrality classes in d + Au collisions. The 
data sets are offset as indicated for clarity, On the top plot is 
the data for and the bottom-plot is the data for 7r° and 
hadron triggered data. The lines shown are a simultaneous 
fit to all the data with Eq. l42l to extract a common fractional 
increase in ^sin 2 (<^jj)^ with T A (b). Bars are statistical errors. 



two increase in CY peripheral to central d + Au colli- 
sions which is much larger than observed in the data. A 
later recombination model by the same authors 23J post- 
dieted only a 30% increase in CY for associated particles 
at px = 2 GeV/c, which is comparable or perhaps slightly 
larger than is observed in the data. 



V. CONCLUSIONS 

We have measured several properties of jet fragmen- 
tation and dijet correlations using two-particle correla- 
tions with three different particle combinations: h ± — h , 
7T° — h , and n — h . From the correlation functions 
we have extracted the widths of the near- and far-angle 
correlations as a function of the momentum of the two 
hadrons, PT.trig and PT.assoc- These widths decrease as a 
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FIG. 41: (Color online) Centrality dependent near (left) and 
far (right) CY(pr) for ti ± — h correlations (upper) and ty° — 
h correlations (lower). The filled circles are from 0-20% 
centrality, the open circles are from 20-40% centrality, the 
filled squares are from 40-88% centrality in d + Au collisions, 
while the open squares are from p + p collisions. Bars are 
statistical errors. 



function of both the trigger and associated particle' s mo- 
menta. From the near-angle widths we calculate v(ir)> 
the RMS of the transverse momentum of fragmented 
hadrons with respect to the hard parton. The value of 
yftjjr) saturates at 0.64 ±0.02 (stat) ±0.04 (sys) GeV/c 
for px, assoc > 2 GeV/c and is consistent with being in- 
dependent of pT.trig and trigger species. The \J (jf,) is 
similar for d + Au and p + p collisions consistent with 
the fragmentation process not being affected by the pres- 
ence of the cold nuclear medium. 

We have also compared the measured x e distributions 
in d + Au collisions to the baseline distributions from 
p + p collisions. The xe distributions extracted from the 
far-angle correlations provide information on the frag- 
mentation of a back-to-back parton triggered on a high- 
Pt hadron in the opposite hemisphere. The measured 
dNf ar /dxE distributions in d + Au are approximately 
independent of pr.trig, i.e. they scale. We have quanti- 
fied the level of scaling violation by extracting the slope 
/3 = d{dN/dxE)/dpT, trig for different ranges of xe- The 
slopes are consistent with zero for d + Au collisions, i.e. 
there is no significant scaling violation. Point-by-point 
the scaling-violation slopes for p + p collisions are not 



FIG. 42: (Color online) Centrality dependence to the ra- 
tio of near (left) and far (right) CY(px) for d + Au to p + p. 

— h correlations are the upper panels and ty° — h correla- 
tions are the lower panels. The shaded band is the systematic 
error due to the normalization of the yields. The filled circles 
are from 0-20% centrality, the open circles are from 20-40% 
centrality, and the filled squares are from 40-88% centrality 
in d+Au collisions. Bars are statistical errors. The shaded 
band represents the normalization error on the spectra. 



significantly different than the d + Au data. This sug- 
gests that if there is any additional gluon radiation in d + 
Au reactions due to multiple scattering, then this has lit- 
tle observable influence on the fragmentation of the hard 
parton. 

We observe no centrality dependence of the conditional 
yield in d + Au and these yields are very similar to 
those from p + p collisions. The recombination model of 
Ref. [23I postdicted a 30% increase in conditional yield 
between d + Au and p+p, which is perhaps slightly larger 
than is observed in the data. 

We have extracted the dijet acoplanarity (sm 2 (cj)jj)) 
from the widths of the back-to-back correlations in d 
+ Au and p + p collisions. In collisions involving 
nuclei, multiple interactions within the nucleus would 
tend to increase the parton transverse momentum which 
would be observable as a larger dijet acoplanarity, i.e., 
the back-to-back distribution of jets should broaden. 
However, in d + Au collisions the extracted values of 
(sin ((pjj)) are very similar to those observed in p + p 
collisions. Indeed, the quadrature difference (A) between 
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(sin 2 (0jj )) in d + Au and p + p is consistent with zero, 
A(sin 2 (0 ij )) = +0.004 ± 0.008(stat)± O.OOO(sys). The 
extracted (sin (<f>jj)) is also measured to be independent 
of the nuclear thickness function Ta(6), which is in con- 
trast to the strongA-dependence of kx observed at lower 
beam energies 0,11, H • 

We have compared the centrality dependence of the 
extracted (sin 2 (4>jjY) with the multiple scattering model 
of Qiu and Vitev jl4| . This model reproduces the mea- 
sured Cronin effect of single-particle spectra at RHIC0] 
3] 4] and predicts a finite increase of kx data with 
nuclear thickness function. When converted to a frac- 
tional increase, the prediction is at a level that is within 
the experimental uncertainty of the current data. Hence 
our present data on (sin 2 (</>jj)) are not inconsistent with 
the level of multiple scattering deduced from the single- 
particle Cronin effect. 

Taken together, we observe no change in fragmentation 
and no indication of the effects of multiple-scattering, i.e. 
the jet-structures are very similar in d+Au and p+p col- 
lisions at RHIC energies. Our measurements also provide 
a critical baseline for jet measurements in Au+Au colli- 
sions at RHIC. 
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